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ABSTRACT 

Numerical models of the tidal disruption of the Sagittarius (Sgr) dwarf galaxy have recently been 
developed that for the first time simultaneously satisfy most observational constraints on the angular 
position, distance, and radial velocity trends of both leading and trailing tidal streams emanating from 
the dwarf. We use these dynamical models in combination with extant 3-D position and velocity data 
for Galactic globular clusters and dSph galaxies to identify those Milky Way satellites that are likely to 
have originally formed in the gravitational potential well of the Sgr dwarf, and have been stripped from 
Sgr during its extended interaction with the Milky Way. We conclude that the globular clusters Arp 
2, M 54, NGC 5634, Terzan 8, and Whiting 1 are almost certainly associated with the Sgr dwarf, and 
that Berkeley 29, NGC 5053, Pal 12, and Terzan 7 are likely to be as well (albeit at lower confidence). 
The initial Sgr system therefore may have contained 5-9 globular clusters, corresponding to a specific 
frequency Sn = 5 — 9 for an initial Sgr luminosity My = —15.0. Our result is consistent with the 8±2 
genuine Sgr globular clusters expected on the basis of statistical modeling of the Galactic globular 
cluster distribution and the corresponding false-association rate due to chance alignments with the 
Sgr streams. The globular clusters identified as most likely to be associated with Sgr are consistent 
with previous reconstructions of the Sgr age-metallicity relation, and show no evidence for a second- 
parameter effect shaping their horizontal branch morphologies. We find no statistically significant 
evidence to suggest that any of the recently discovered population of ultra-faint dwarf galaxies are 
associated with the Sgr tidal streams, but are unable to rule out this possibility conclusively for all 
systems. 

Subject headings: Sagittarius Dwarf - Milky Way: globular clusters - Milky Way: structure - Local 
Group 
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1. INTRODUCTION 

Globular clusters tend to form in all sufficiently mas- 
sive galactic systems during episodes of major star forma- 
tion, their relatively simple stellar populations providing 
important constraints on starburst and galaxy formation 
models (see, e.g., Brodie & Strader 2006). While no 



globular clusters are found in low mass (M\, 
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dwarf galaxies, the more massive dwarfs such as the 
LMC (M v = -18.5), SMC (My = -17.1) and Fornax 
(My = —13.1) are observed to contain nineteen, eight, 
and five globular clusters respectively (see, e.g., Forbes 
et al. 2000). It is generally accepted that the Sagittar- 
ius (Sgr) dwarf spheroidal galaxy {My = —13.64; Law & 
Majewski 2010 [hereafter LM10]) contains at least four 
clusters (M 54, Arp 2, Terzan 7, Terzan 8; Da Costa 
& Armandroff 1995) within its main body. Since Sgr 
is known to have experienced significant tidal mass loss, 
however, it is likely that additional globular clusters may 
have been stripped from Sgr during its prolonged interac- 
tion with the Milky Way and now lie scattered through- 
out the Galactic halo. 

In the anticipation that such clusters are likely dis- 
tributed across the sky, various authors (e.g., Lynden- 
Bell & Lynden-Bell 1995; Irwin 1999; Bellazzini et al. 
2002; Newberg et al. 2003; Majewski et al. 2004; Car- 
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raro et al. 2007) have suggested a number of additional 
globular clusters (e.g., Pal 2, Pal 12, NGC 5634, NGC 
2419, Whiting 1, etc.) that may have formed in the Sgr 
dSph (see Table Q] for a summary). Perhaps the most 
comprehensive efforts to conduct a systematic census to 
date have been undertaken by Palma et al. (2002) and 
Bellazzini et al. (2003a), who demonstrated that there 
is a strong correlation between the positions and radial 
velocities of outer halo (Galactocentric radius tqc > 10 
kpc) globular clusters and the orbit of Sgr. Bellazzini et 
al. (2003a) in particular presented a statistical assess- 
ment of the significance of the total number of Galac- 
tic globular clusters aligned with the Sgr orbital plane, 
concluding that there was a < 2% probability that the 
observed grouping of clusters about the Sgr plane arises 
by chance. While this approach was certainly effective, 
it was based on the projected orbit of Sgr (as derived by 
Ibata & Lewis 1998) and therefore did not adopt an op- 
timal membership criterion because tidal debris does not 
precisely follow the orbit of the parent object: Trailing 
tidal arms lie outside the orbital path of the parent while 
leading tidal arms lie interior to the orbit and, in the case 
of Sgr, wrap up around the Galactic Center (see discus- 
sion by Johnston et al. 1995, 1999; Choi et al. 2007; 
Eyre & Binney 2009). Moreover, Bellazini et al. (2003a) 
concentrated on the past orbit of Sgr, an approach that 
misses potential connections to the leading arm tidal de- 
bris. 

In this contribution we wish to update and extend 
these previous analyses, and ask whether each individual 
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Milky Way satellite matches the angular position, dis- 
tance, and radial velocity of Sgr tidal debris sufficiently 
well that it is statistically likely to be physically associ- 
ated with the stream and might therefore have originated 
in the Sgr dwarf. Such an undertaking is greatly aided 
by the significantly improved picture of the Milky Way 
— Sgr system that has been provided in recent years by 
the Two Micron All-Sky Survey (2MASS) and Sloan Dig- 
ital Sky Survey (SDSS) which have mapped the stellar 
streams from the Sgr dwarf wrapping fully 360° across 
the sky (Majewski et al. 2003; Bclokurov et al. 2006a; 
Yanny et al. 2009). These surveys have also revealed a 
wealth of additional substructure in the Galactic halo in 
the form of both globular clusters and a population of 
dwarf galaxies with masses ~ 1Q 7 Mq within their cen- 
tral 300 pc (Strigari et al. 2008) but ultra-faint luminosi- 
ties comparable to those of globular clusters. The origin 
of these ultra-faint dwarfs is unknown, and we consider 
whether any of them may be dynamically associated with 
the Sgr dSph. 

Our effort is also aided by the recent construction of 
a numerical model of the Sgr stream that reproduces 
the majority of the observed characteristics of the stel- 
lar tidal streams with extremely high fidelity (Law et al. 
2009; LM10). The comprehensive view of both leading 
and trailing arm material (which overlap each other along 
the line of sight for large swathes of the sky) provided by 
this model allows us to view potential Sgr members in 
the context of the full debris system, and assess possi- 
ble associations with leading and trailing arms based not 
only on spatial distributions, but on dynamical proper- 
ties such as radial velocities and proper motions as well. 

This paper is organized as follows. In Sj2]we describe 
the general characteristics of the Sgr tidal streams and 
the iV-body model that best reproduces its observed 
properties. In Sj3] we describe our sample of Galac- 
tic satellites, incorporating globular clusters, ultra-faint 
dwarfs, and two open clusters for which association with 
Sgr has previously been claimed in the literature. In £|4] 
we describe our numerical technique for quantifying the 
association of a Galactic satellite with the Sgr stream, 
applying this to individual stellar clusters and ultra-faint 
dwarfs in fj5]and ^respectively. We discuss the result- 
ing implications for the star cluster budget of Sgr and 
its contributions to the Galactic halo in SJ71 summarizing 
our conclusions in EjSj 

We adopt the heliocentric Sgr coordinate system 
(A©,_B Q ) defined by Majewski et al. (2003) in which 
the longitudinal coordinate A© = 0° in the direction of 
Sgr and increases along trailing tidal debris, and B & is 
positive towards the orbital pole (I, b) po i c rs (274°, —14°) 
(see discussion by LM10). All radial velocities are given 
in the Galactic Standard of Rest (GSR) frame, with re- 
spect to which the Sun has a peculiar motion (U, V, W) 
= (9, 12 + 220, 7) km s _1 and is located R Q = 8 kpc from 
the Galactic Center. 

2. THE SAGITTARIUS DISRUPTION MODEL 

We elect to compare the Galactic satellite population 
to the LM10 iV-body modeQ of the Sgr streams. As dis- 
cussed at length in LM10, this model was constrained to 



4 This model is available online 
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match the abundant observational data on the angular 
positions (e.g., Majewski et al. 2003; Belokurov et al. 
2006a; Yanny et al. 2009), radial velocities (e.g., Law et 
al. 2004, 2005; Majewski et al. 2004; Monaco et al. 2007; 
Yanny et al. 2009), and chemical abundances (Chou et 
al. 2007, 2010; Monaco et al. 2007; Yanny et al. 2009) of 
stars in the Sgr tidal streams which have recently been 
provided by 2MASS and SDSS. This model was fully con- 
strained using only angular positions and radial veloci- 
ties for Sgr stream stars; photometric parallaxes were not 
used as a constraint because of their large observational 
uncertainties and systematic effects due to variations in 
metallicity and age along the tidal arms (e.g., Chou et 
al. 2007, 2010). As detailed in LM10, however, the N- 
body model matches distance estimates for stars in the 
Sgr streams to within observational uncertainty. 

We fit to the LM10 TV-body model instead of the raw 
observational data for four reasons: (1) The positional 
data for the Sgr M-giant stream (Majewski et al. 2003) 
contains an intrinsic distance scatter (estimated to be 
17% by Majewski et al. 2003) from the intrinsic width of 
the color-magnitude relation, as well as poorly known 
systematic effects in distance because of variations in 
the metallicity distribution function of Sgr stream stars. 
Similar systematic uncertainties also affect the SDSS ob- 
servations, making certain sections of the tidal arms dif- 
ficult to convert to distance directly. The LM10 model 
incorporates these metallicity effects and is fully consis- 
tent with the M-giant color-magnitude relation and ap- 
parent magnitude trend. (2) The available velocity data 
are still incomplete, with whole sections of the tidal arms 
poorly sampled. The LM10 model can fill in these gaps 
in our knowledge. (3) There is a limit to the length of the 
Sgr tidal tails that can be traced with individual stellar 
populations such as M-giants (Majewski et al. 2003) or 
blue horizontal branch stars (Yanny et al. 2009). Ad- 
ditional tidal debris from older stellar generations may 
exist and simply be poorly traced by current observa- 
tional data. With the LM10 model, the observed tails 
may be extrapolated to include such hypothetical older 
tidal debris. (4) The LM10 model can provide proper 
motion information at all positions along the Sgr tidal 
streams. These proper motions are well-constrained by 
observations in the other four dimensions of phase space, 
but are extremely difficulty to measure directly. 

In Figure Q] (colored points) we plot a cross-sectional 
view of the LM10 model of the Sgr system in Galacto- 
ccntric Cartesian coordinates. The model satellite has 
been orbiting in the Galactic potential for ~ 8 Gyr, and 
has produced multiple epochs of tidal debris correspond- 
ing to each of its pericentric approaches to the Milky 
Way. We adopt the LM10 scheme in which tidal debris 
is color-coded according to the pair of orbits on which 
it became unbound from the Sgr core: Black points are 
currently bound to Sgr, orange points became unbound 
in the last 1.3 Gyr, magenta points between 1.3 and 3.2 
Gyr ago, cyan points between 3.2 and 5.0 Gyr ago, and 
green points between 5.0 and 6.9 Gyr ago (see Fig. 7 of 
LM10). 

Since the Sgr streams frequently overlap themselves 
in angular position, distance, and radial velocity, it is 
important to distinguish the individual wraps of these 
streams from each other when evaluating possible asso- 
ciations with stellar subsystems in the Galactic halo. We 
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TABLE 1 

Candidate Sgr Stream Globular Clusters Proposed in the Literature 81 



ID 



Name 



DA95 



LL95 



199 



D01 



P02 



B03 



M04 



C07 



LMlOb 



■15 
52 

21 
42 
1 
13 

18 
22 
2(5 
27 
30 
10 

7 

31 

53 

51 
19 

11 
46 
3 



Arp 2 
M 2 
M 5 
M 53 
M 54 
NGC 288 
NGC 2419 
NGC 4147 
NGC 5053 
NGC 5466 
NGC 5634 
NGC 5824 
NGC 6426 
NGC 6356 
Pal 2 
Pal 5 
Pal 12 
Pal 13 
Rup 106 
Terzan 3 
Terzan 7 
Terzan 8 
Whiting 1 



YES 



YES 



YES 



YES 



YES 

YES 
YES 



NO 



YES 


YES 




YES 


r ES YES 


YES 


'ES NO 


YES 



YES 


YES 
YES 


YES 


YES 






YES 






YES 


YES 


YES 




YES 




NO 


YES 




NO 


YES 




YES 


YES 
YES 
YES 
YES 




YES 








YES 


YES 


NO 


YES 




YES 


YES 


YES 


NO 


YES 
YES 
YES 




YES 


YES 


YES 


YES 


YES 


YES 



YES 



YES 
NO 
NO 
NO 
YES 
NO 
NO 
NO b 
YES 
NO 
YES 
NO 
NO 
NO 
NO b 
NO 
YES 
NO 
NO 
NO 
YES 
YES 
YES 



a Columns represent Da Costa & Armandroff (1995; DA95), Lynden-Bcll & Lynden-Bcll (1995; L95), Irwin et al. (1999; 199), Dinescu et al. (2001; 
D01), Palma et al. (2002; P02), Bellazzini et al. (2003a; B03), Majewski et al. (2004; M04), Carraro et al. (2007; C07), this contribution (LMlOb). 
'YES' indicates possible candidate, 'NO' indicates disfavored by indicated study. 

As discussed in ^\ NGC 4147 and Pal 2 may plausibly be associated with the Sgr stream, but at low statistical confidence. 



therefore follow LM10 in adopting a shorthand notation 
describing specific wraps of the streams (see also Fig. 8 
ofLMlO): 

LI:: Primary wrap of the leading arm of Sgr (i.e., 0° — 
360° angular separation from the dwarf). Traced 
by orange/magenta points near the Sgr dwarf, ma- 
genta/cyan/green points at larger angular separa- 
tions. 

L1Y:: Young subset of the LI wrap, incorporating only 
orange/magenta points (i.e., debris lost during the 
last 3 Gyr). 

L2:: Secondary wrap of the leading arm of Sgr (i.e., 
360° — 720° angular separation from the dwarf). 
Traced almost entirely by cyan/green points. 

Tl:: Primary wrap of the trailing arm of Sgr (i.e., 0° — 
360° angular separation from the dwarf). Traced 
by orange/magenta points near the Sgr dwarf, ma- 
genta/cyan/green points at larger angular separa- 
tions. 

T1Y:: Young subset of the Tl wrap, incorporating only 
orange/magenta points (i.e., debris lost during the 
last 3 Gyr). 

T2:: Secondary wrap of the trailing arm of Sgr (i.e., 
360° — 720° angular separation from the dwarf). 
Traced almost entirely by cyan/green points. 

The L1Y and T1Y stream segments described above 
correspond most closely to those phases of tidal debris 
< 3 Gyr in age for which there is the strongest obser- 
vational evidence from the 2MASS and SDSS wide-field 



surveys. The true length of the debris streams is un- 
known however, and it is possible that they may extend 
substantially further than has been conclusively observed 
to date. The L2/T2 streams (and the cyan/green sec- 
tions of Ll/Tl) represent model predictions for where 
such putative older tidal debris may lie, should it be 
present in the Galactic halo. 

3. SELECTING A SATELLITE POPULATION 

As illustrated in Figure [TJ the most recent (i.e., < 3 
Gyr in age) Sgr tidal debris lies at distances 10 < d < 60 
kpc from the Sun, and greater than 8 kpc from the Galac- 
tic Center. The LM10 simulations suggest, however, that 
putative L2/T2 tidal debris may exist as much as 100 
kpc away. We therefore construct our comparison sam- 
ple of Galactic globular clusters by extracting all those 
with Galactocentric radii 8 < tqc < 100 kpc from the 
Harris (1996, with Feb. 2003 update) catalog; this con- 
servatively includes all well-known globular clusters that 
might conceivably be associated with Sgr. We note that 
this rules out checking for possible angular alignments 
with much more distant clusters that may have been lost 
extremely at extremely early times; we discuss this point 
further in §5.191 Motivated by claims for association with 
Sgr tidal debris by Carraro et al. (2007, 2009), we also 
include in our sample the star clusters Whiting 1, Berke- 
ley 29 (Be 29), and Saurer 1 (Sa 1). While Be 29 and Sa 
1 are commonly regarded as Galactic open clusters, their 
physical properties overlap with the lower-mass globular 
clusters in our sample, and we consider that they may 
be sparse halo globular clusters which happen to lie near 
the Galactic disk. 

Recent years have witnessed the discovery in the Sloan 
Digital Sky Survey (SDSS) of a multitude of "ultra-faint" 
Galactic satellites that have significantly lower luminos- 
ity than the classic dwarf spheroidal galaxy population, 
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Fig. 1. — Diagram showing the LM10 model of the Sgr system (colored points) in Galactocentric Cartesian coordinates (XqCi-^Gc)- 
The color of individual debris particles represents the time at which the particle became unbound from the Sgr dwarf. The X indicates the 
location of the Galactic Center, while the circled X indicates the location of the Sun. The horizontal dashed line represents the Galactic 
disk plane, the length of the line is chosen to indicate a diameter 30 kpc. Overplotted against the model are ID numbers corresponding to 
the locations in the Xqq — Zqc plane of Milky Way stellar satellites (ID numbers are matched with the corresponding satellite names in 
Tabled- The locations of some ID numbers have been adjusted slightly to minimize confusion from overlapping labels. 



but half-light radii too large to be traditional globular 
clusters (see, e.g., discussion by Gilmore et al. 2007). 
The origin of this low luminosity population is uncer- 
tain, and some have been suggested to be stellar over- 
densities previously associated at some time in the past 
with massive dSphs such as Sgr. We therefore include 
in our analysis the ultra- faint systems Bootes I (Boo I), 
Bootes II (Boo II), Coma Berenices (Coma Ber), Pisces 
I, Segue 1, Segue 2, Segue 3, Ursa Major II (UMa II), and 
Willman I. The Canes Venatici I, Canes Venatici II, Her- 
cules, Leo IV, Leo V, Leo T, and Ursa Major I systems 
along with the classical dSph population (e.g., Carina, 
Draco, Fornax, etc.) are excluded from formal consider- 
ation since they lie at large distances (typically greater 
than 100 kpc) at which there is no Sgr tidal debris in the 
LM10 model, although we discuss some of them briefly 



ingfUOl 

The final sample of 64 satellites is tabulated in Table 
[2] (sorted in order of increasing right ascension); 59 of 
these satellites have both distance and radial velocity 
measurements published in the literature. Of these 59 
satellites, 51 are stellar clusters and eight are ultrafaint 
dwarfs. Proper motions^ have been measured for 24 of 
these satellites, and have been drawn from the Harris 
(1996) catalog, Dinescu et al. (1999, 2001, 2003), Siegel 
et al. (2001), Palma et al. (2002), and Cassetti-Dinescu 
et al. (2007). 

5 Because the apparent proper motion of the Sgr stream along 
the direction of Galactic longitude I varies extremely rapidly for 
tidal debris passing near the poles of the Galactic (I, b) coordinate 
system we opt to present proper motions in the (a, 8) coordinate 
frame. 
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TABLE 2 
Properties of the Milky Way Satellites 



ID 


Name 


R.A. 


Dccl. 


I 


b 


l-^Sgrl 


A 


B Q a 


Distance 


Vgsr 


fia,cosS 


f-s 






(J2000.0) 


(J2000.0) 


(°) 


(°) 


(kpc) 


(°) 


(°) 


(kpc) 


(km s _1 ) 


(mas yr _1 ) 


(mas yr _1 ) 


Star Clusters 


1 


NGC 288 


00:52:48 


-26:35:24 


152 


-89 


1.1 


76 


13 


8.8 


-52.5 + 0.4 


4.68 + 0.18 


-5.37 + 0.52 


2 


NGC 362 


01:03:14 


-70:50:54 


302 


-46 


5.59 


46 


50 


8.5 


85 + 0.5 


5.08 + 1.15 


-2.49 + 2.17 


3 


Whiting 1 


02:02:57 


-03:15:10 


162 


-61 


0.215 


103 


1 


29.4 


-105 + 1.8 






1 


NGC 1261 


03:12:15 


-55:13:01 


271 


-52 


11.9 


79 


51 


16.4 


-79.7 + 4.6 






5 


Pal 1 


03:33:23 


+79:34:50 


130 


19 


9.8 


180 


-55 


10.9 


81.9 + 3.3 






6 


Eridanus 


04:24:45 


-21:11:13 


218 


-41 


50.1 


125 


35 


90.2 


-141 ±2.1 






7 


Pal 2 


04:46:06 


+31:22:51 


171 


-9 


5.95 


154 


-11 


27.6 


-105 + 57 






8 


NGC 1851 


05:14:06 


-40:02:50 


245 


-35 


9.13 


121 


56 


12.1 


142 ± 0.6 


1.28 + 0.68 


2.39 ± 0.65 


9 


M 79 (NGC 1904) 


05:24:11 


-24:31:27 


227 


-29 


8.1 


138 


44 


12.9 


48.8 ± 0.4 


2.12 ± 0.64 


-0.02 ± 0.64 


10 


NGC 2298 


06:48:59 


-36:00:19 


246 


-16 


8.61 


157 


63 


10.7 


-59.7 + 1.2 


4.05 ± 1 


-1.72 + 0.98 


11 


Berkeley 29 


06:53:04 


+16:55:39 


198 


8 


1.79 


177 


12 


13.2 


-53.7 + 3.6 






12 


Saurer 1 


07:20:54 


+01:48:00 


215 


7 


5.25 


182 


28 


13.2 


-32.8 + 3.6 






13 


NGC 2419 


07:38:09 


+38:52:55 


180 


25 


13.6 


190 


-9 


84.2 


-26.5 + 0.8 






14 


Pyxis 


09:07:58 


-37:13:17 


261 


7 


35.4 


224 


66 


39.7 


-194 + 1.9 






15 


NGC 2808 


09:12:03 


-64:51:47 


282 


-11 


8.6 


332 


82 


9.6 


-128 + 2.4 


0.58 + 0.45 


2.06 ± 0.46 


16 


Pal 3 


10:05:31 


+00:04:17 


240 


42 


40.6 


228 


27 


92.7 


-65.1 + 8.4 


0.33 + 0.23 


0.3 + 0.31 


17 


NGC 3201 


10:17:37 


-46:24:40 


277 


9 


3.73 


265 


68 


5 


269 ± 0.2 


5.28 + 0.32 


-0.98 + 0.33 


18 


NGC 4147 


12:10:06 


+18:32:31 


253 


77 


1.39 


251 


-1 


19.3 


140 ± 0.7 


-1.63 + 1.44 


-1.86 + 3.21 


19 


Rup 106 


12:38:40 


-51:09:01 


301 


12 


16.1 


304 


53 


21.2 


-233 ± 3 






20 


M 68 (NGC 4590) 


12:39:28 


-26:44:34 


300 


36 


4.92 


281 


35 


10.2 


-250 + 0.4 


-3.76 + 0.66 


1.79 + 0.62 


21 


M 53 (NGC 5024) 


13:12:55 


+18:10:09 


333 


80 


3.4 


265 


-8 


17.8 


-89.5 + 4.1 


0.5 + 1 


-0.1 + 1 


22 


NGC 5053 


13:16:27 


+17:41:53 


336 


79 


3.21 


266 


-8 


16.4 


34.1 + 0.4 






23 


M 3 (NGC 5272) 


13:42:11 


+28:22:32 


42 


79 


4.5 


265 


-20 


10.4 


-109 + 0.2 


-2.84 + 2.01 


-2.45 + 0.62 


21 


NGC 5286 


13:46:27 


-51:22:24 


312 


11 


6.94 


315 


45 


11 


-106 + 1.5 






25 


AM 4 


13:55:50 


-27:10:22 


320 


34 


12 


298 


26 


29.9 








26 


NGC 5466 


14:05:27 


+28:32:04 


42 


74 


7.16 


270 


-23 


15.9 


160 + 0.3 


-4.46 ± 1.32 


0.85 ± 0.84 


27 


NGC 5634 


14:29:37 


-05:58:35 


342 


49 


0.531 


293 


3 


25.2 


-80.5 + 6.6 






28 


NGC 5694 


14:39:37 


-26:32:18 


331 


30 


10.7 


306 


20 


34.7 


-231 + 1.3 






29 


IC 4499 


15:00:19 


-82:12:49 


307 


-20 


15 


2.89 


57 


18.9 








30 


NGC 5824 


15:03:59 


-33:04:04 


333 


22 


11.2 


315 


22 


32 


-117 + 1.5 






31 


Pal 5 


15:16:05 


+00:06:41 


1 


46 


4.06 


300 


-8 


23.2 


-44.3 + 0.2 


-2.29 + 0.67 


-1.65 + 0.81 


32 


BH 176 


15:39:07 


-50:03:02 


328 


4 


7.58 


332 


33 


15.6 








33 


Pal 14 


16:11:05 


+14:57:29 


29 


42 


35 


305 


-27 


73.9 


170 + 2.2 






34 


NGC 6101 


16:25:49 


-72:12:06 


318 


-16 


10.4 


355 


47 


15.3 


216 + 1.7 






35 


M 13 (NGC 6205) 


16:41:42 


+36:27:37 


59 


41 


6.72 


297 


-49 


7.7 


-87.2 + 0.3 


-0.9 + 0.71 


5.5 + 0.89 


36 


NGC 6229 


16:46:59 


+47:31:40 


74 


40 


26.6 


285 


-58 


30.4 


22 ± 7.6 






37 


Pal 15 


17:00:02 


+00:32:31 


20 


25 


16.2 


324 


-20 


44.6 


151 + 1.1 






38 


M 92 (NGC 6341) 


17:17:07 


+43:08:11 


68 


35 


7.9 


298 


-59 


8.2 


63.4 + 0.1 


-3.75 + 1.72 


0.39 + 1.68 


39 


IC 1257 


17:27:09 


-07:05:35 


17 


15 


7.61 


334 


-16 


25 


-66.3 + 2.1 






40 


NGC 6426 


17:44:55 


+03:10:13 


28 


16 


10.2 


334 


-27 


20.7 


-47.5 ± 23 






41 


ESO 280-SC06 


18:09:06 


-46:25:23 


347 


-13 


6.17 


355 


19 


21.7 








42 


M 54 (NGC 6715) 


18:55:03 


-30:28:42 


6 


-14 


0.199 





1 


28.0 


171 + 0.5 






43 


M 56 (NGC 6779) 


19:16:36 


+30:11:05 


63 


8 


9.56 


352 


-59 


10.1 


73.3 + 0.8 


0.3 + 1 


1.4+1 


11 


Terzan 7 


19:17:44 


-34:39:27 


3 


-20 


1.09 


5 


5 


23.2 


185 + 4 






45 


Arp 2 


19:28:44 


-30:21:14 


9 


-21 


0.71 


7 





28.6 


153 ± 10 






46 


Terzan 8 


19:41:45 


-34:00:01 


6 


-25 


0.812 


10 


4 


26 


156 + 8 






47 


M 75 (NGC 6864) 


20:06:05 


-21:55:17 


20 


-26 


3.97 


15 


-9 


20.7 


-112 + 3.6 






48 


NGC 6934 


20:34:12 


+07:24:15 


52 


-19 


10.5 


23 


-38 


15.7 


-235 + 1.6 


1.2 + 1 


-5.1 + 1 


49 


M 72 (NGC 6981) 


20:53:28 


-12:32:13 


35 


-33 


6.01 


27 


-18 


17 


-230 + 3.7 






50 


NGC 7006 


21:01:30 


+16:11:15 


64 


-19 


30.7 


33 


-46 


41.5 


-186 + 0.4 


-0.96 ± 0.35 


-1.14 + 0.4 


51 


M 15 (NGC 7078) 


21:29:58 


+ 12:10:01 


65 


-27 


7.6 


42 


-41 


10.3 


80 + 0.2 


-1.28 + 1.34 


-7.05 ± 3.09 


52 


M 2 (NGC 7089) 


21:33:29 


+00:49:23 


55 


-35 


6.53 


40 


-29 


11.5 


151 + 2 


6.03 ± 0.99 


-5.19 + 1.34 


53 


Pal 12 


21:46:39 


-21:15:03 


31 


-48 


3.24 


39 


-7 


19.1 


107+1.5 


-1.2 + 0.3 


-4.21 + 0.29 


51 


Pal 13 


23:06:44 


+12:46:19 


87 


-43 


15.1 


70 


-34 


25.8 


190 + 0.5 


2.3 + 0.26 


0.27 + 0.25 


55 


NGC 7492 


23:08:27 


-15:36:41 


53 


-63 


4.06 


59 


-7 


25.8 


-128 + 4.3 






Ultrafaint Dwarfs 


56 


Segue 2 


02:19:16 


+20:10:31 


149 


-38 


11 


119 


-17 


35 


43.2 + 2.5 






57 


UMall 


08:51:30 


+63:07:48 


152 


37 


18.3 


202 


-33 


32 


-33.4+1.9 






58 


Segue 1 


10:07:03 


+16:04:25 


220 


50 


3.5 


225 


11 


23 


111 + 1.3 






59 


Willman I 


10:49:22 


+51:03:04 


159 


57 


16.9 


223 


-25 


38 


35.4 + 2.5 






60 


Coma Bcr 


12:26:59 


+23:54:15 


242 


84 


7.03 


252 


-8 


44 


81.8 + 0.9 






61 


Bootes II 


13:58:05 


+12:52:00 


354 


69 


8.09 


277 


-9 


46 


-116 + 5.2 






62 


Bootes I 


14:00:06 


+14:30:00 


358 


70 


12.4 


277 


-11 


62 


103 + 3.4 






63 


Segue 3 


21:21:31 


+19:07:02 


69 


-21 


12.8 


41 


-48 


16 








64 


Pisces I 


23:40:00 


+00:18:00 


88 


-58 


28.2 


73 


-19 


85 


45.4 + 4 
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4. QUANTIFYING THE ASSOCIATION OF SATELLITES 

While the best evaluation of the association of any indi- 
vidual satellite can be made by considering the properties 
of the satellite in comparison to those of the Sgr stream 
in detail, a preliminary cut of reasonable stream candi- 
dates can be made via simple statistical arguments. In 
this section, we describe our construction of a correlation 
statistic that can be used to pre-select those Milky Way 
satellites that match the Sgr streams sufficiently well to 
warrant close inspection in $5] and [SI 

4.1. Defining an Association Statistic 

We quantify the offset of each satellite listed in Tabled 
from the LM10 model of the Sgr stream as the quadrature 
sum of the angular separation, difference in heliocentric 
distance, and difference in radial velocity between the 
satellite and the stream: 



X 



(B e - P Q , Sgr ) 2 (d - d S gr) 2 , (V - W S gr) 



(1) 



M.Sgr 



v,Sgr 



where Bq , d, and v are the latitude (in Sgr plane coordi- 
nates), the heliocentric distance, and the radial velocity 
of the satellite respectively. B©,sgD 0\b Sgr , efegr, Cd,Sgr, 
Uggr, and tr V] g gr respectively represent the mean and la 
width of the latitude, distance, and radial velocity of the 
LM10 stream, calculated within 5° of orbital longitude 
A© of the satellite (i.e., taking A Q as the independent 
variable). We assume a typical observational uncertainty 
Pd/d = 10%, and a v = 3 km s _1 for each of the satellites. 
X 2 therefore quantifies the offset between a satellite and 
the Sgr stream, weighted by the quadrature sum of the 
observational uncertainty and the variance of Sgr stream 
stars about the mean at the corresponding location. \ 2 
is calculated with respect to each of the arms of the Sgr 
stream listed above in 82] (i.e., x L i> Xliy> Xl,2> etc 0- 

4.1.1. Application to the Globular Clusters 

Taken on its own, \ 2 is a marginally useful statistic 
that quantifies how closely a given satellite matches the 
properties of the Sgr stream. A more interesting ques- 
tion, however, is the following: How significant is the 
apparent association compared to that expected for a 
population of satellites randomly distributed through- 
out the Galactic halo? To assess this, we perform a 
series of Monte Carlo simulations, randomly populating 
the Galactic halo with 10 4 artificial satellites between 
Galactocentric radii 8 kpc < tqc < 100 kpc. We as- 
sume a spherically symmetric distribution with a radial 
density profile proportional to Tqq (Bellazzini et al. 
2003a), which mimics the Galactic halo globular cluster 
population in this range of Galactocentric radii. These 
physical locations are transformed to the heliocentric 
(Aq, BQ,d) coordinate frame, and paired with a radial 
velocity drawn at random from a Gaussian distribution 
with mean («gsr) = — 38 km s _1 and standard deviation 
tr^ GSR = 175 km s" 1 (see Bellazzini et al. 2003a)Q 

Calculating % 2 for each of these 10 4 artificial globular 
clusters, we determine the fraction P(x) that by chance 
match some arm of the Sgr stream to an accuracy of \ 
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6 Adopting (vgsr) = km s 
ference to our calculations. 



1 instead makes a negligible dif- 



FlG. 2. — Probability P(x) for an artificial halo globular cluster 
to match the Sgr stream with an accuracy of \ or better. The 
solid line is the relation obtained when considering only matches 
to the most recent 3 Gyr of tidal debris (i.e., L1Y and T1Y wraps), 
the dashed line is the relation obtained when considering matches 
to the most recent 7 Gyr of tidal debris (i.e., LI, Tl, L2, and T2 
wraps). 

or better. Since the likelihood of matching some seg- 
ment of the Sgr stream depends strongly on the length 
of the stream (if the streams were arbitrarily long they 
would fill the available parameter space, and any ran- 
domly distributed satellite would match some wrap rea- 
sonably well) we calculate two versions of the P statistic: 

P3:: Fraction of randomly distributed satellites that 
match either the L1Y or T1Y arms (i.e., the last 3 
Gyr of tidal debris) to an accuracy better than %. 

P7:: Fraction of randomly distributed satellites that 
match either the LI, L2, Tl, or T2 arms (i.e., the 
last 7 Gyr of tidal debris) to an accuracy better 
than x- 

If Sgr has only been interacting with the Milky Way 
for the last 3 Gyr (i.e., there is no significant tidal debris 
beyond that conclusively traced to date by 2MASS and 
SDSS) the P3 statistic will be most realistic, whereas if 
there is appreciable tidal debris from earlier epochs not 
traced by current surveys then the P7 statistic will be 
most appropriate. We plot P3 and P7 as functions of 
X in Figure [2] noting that P3 and P7 can also be in- 
terpreted as the probability that some randomly chosen 
artificial globular cluster will match the Sgr streams to 
an accuracy of x or better. Figure [2] illustrates that there 
is a low probability for a randomly chosen artificial clus- 
ter to match the Sgr stream to high accuracy (x < 3), 
while nearly all artificial clusters will match the stream 
to better than \ ~ 30. Unsurprisingly, P7 rises faster 
than P3 as a function of x because there are more oppor- 
tunities for each cluster to match some wrap of the tidal 
streams. 

Using the relations shown in Figure [2] it is possible 
to calculate the probabilities P7.L1, P3.LIY, etc., that a 
randomly chosen artificial cluster would match some seg- 
ment of the Sgr stream as well as (or better than) the 
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Fig. 3. — Number of globular clusters as a function of the proba- 
bility statistic P. A randomly distributed cluster population would 
be uniform from P = — 1 (solid horizontal line) modulo statisti- 
cal fluctuations, while the actual Milky Way distribution exhibits 
a peak at P < 0.15 on account of the presence of clusters genuinely 
associated with the Sgr stream. The red histogram corresponds to 
association probabilities considering all wraps of the LM10 model 
stream (P7), the blue histogram indicates the result when only 
dynamically young sections most closely corresponding to observa- 
tional data are used (P3). 

real cluster matches a given segment of the stream. Val- 
ues of these P are tabulated for each of our 51 globular 
clusters with full angular, distance, and radial velocity 
information in Tables [3] and 01 Note that P3 is defined 
as the minimum of P3.LIY and i^TiY, and Pj as the 
minimum of Pj,li, Ptjji, -FV,ti, and Pj, T2- P3 and P7 
therefore summarize the significance of the association of 
a given cluster with the Sgr streams under two different 
assumptions for the length of the streams. 

In Figure [3] we plot a histogram of P3 and Pr for 
our 51 globular clusters. If these globular clusters were 
distributed randomly throughout the Galactic halo, the 
number N(P) of satellites in each bin would be uni- 
form throughout the range P = — 1 (modulo statistical 
fluctuations). The spike in the number of clusters with 
P < 0.15 reflects the fact that there is a globular clus- 
ter population genuinely associated with the Sgr tidal 
streams; we therefore adopt the criterion P3 < 0.15 or 
P-i < 0.15 for selecting globular clusters as candidates 
for association with Sgr. Based on our Monte Carlo 
simulations of 10 4 randomly distributed artificial clus- 
ters, and 500 realizations of randomly extracted sets of 
51 such clusters, we expect that 7 ± 2 globular clusters 
will meet the selection criteria (either P3 or Pj) based on 
chance alignment with the Sgr streams. In comparison, 
15/14 globular clusters actually meet the P3/P7 criteria 
respectively^] In less than 1% of realizations of a sam- 
ple of 51 randomly drawn artificial clusters do we find 
such a large number with P3 or P7 < 0.15 by chance. 

7 Note that some clusters that are Sgr stream candidates using 
the P3 statistic are not candidates under the Pj statistic because 
by opening the parameter space for potential matches to include 
alignment with the older arms, the Monte Carlo probability nor- 
malization necessarily imposes more stringent constraints on pos- 
sible associations with the young arms. 
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Fig. 4. — As Figure[21 but for an artificial population constructed 
to mimic the distribution of Milky Way ultra-faint satellites. 




Fig. 5. — As Figure[3] but for the ultra- faint Milky Way satellites. 

Subtracting from the maximum number of potential Sgr 
globular clusters (15) the typical number expected to be 
false associations (7±2), we conclude that there are likely 
8 ± 2 globular clusters genuinely associated with the Sgr 
stream. 

4.1.2. Application to the Ultrafaint Satellites 



The P statistic defined in §4.1.11 above is not appro- 
priate for the ultra-faint satellite population discovered 
in the SDSS. Not only do these satellites have a different 
radial number density profile than the globular clusters, 
but they are subject to a strong selection bias in the sense 
that only satellites lying within the SDSS footprint could 
have been discovered. Since the Sgr stream occupies a 
significant fraction of this footprint, the ultra-faint satel- 
lites discovered to date will therefore be predisposed to 
lie relatively close to the stream. 

We therefore develop an alternative formulation of the 
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P statistic appropriate to the ultra-faint satellites. Since 
the radial number density profile of this population is 
poorly constrained by the small number of such satel- 
lites discovered to date, we generate an artificial compar- 
ison population whose Galactocentric radii match those 
of the observed satellites. We assume that this artifi- 
cial population is spherically symmetric about the Milky 
Way, but reject from consideration any satellites that fall 
outside of the SDSS footprint^ Sufficiently many artifi- 
cial satellites are generated that the final population that 
falls within the SDSS footprint is composed of 10 4 mem- 
bers, for each of which a radial velocity is chosen from 
a Gaussian random distribution with mean (i>gsr) = 
km s^ 1 and standard deviation o~ VGSIS = 75 km s _1 that 
fits the observed distribution of wgsr for the real satel- 
lites in question. 

In Figure [4] we plot P as a function of x f° r this ar ~ 
tificially generated comparison population, noting that 
its cumulative distribution rises much more steeply than 
did the cumulative distribution of the artificial globu- 
lar clusters (Figure [2]) . This reflects the fact that since 
the ultra-faint satellites by definition lie within the SDSS 
footprint, they are predisposed to relatively strong ap- 
parent correlations with the Sgr stream. Figure [5] plots 
the histogram of P$ and Pj for the eight ultra- faint satel- 
lites whose radial velocities have been measured. The 
roughly uniform distribution between P = — I indicates 
that the satellites are statistically consistent with being 
unassociated with the Sgr stream. In particular, there 
is no clump at low values of P corresponding to close 
association with the Sgr streams. Indeed, there may be 
weak evidence (at the ~ 2er level) using the P3 statistic 
for an anti- correlation between some of these satellites 
and the Sgr streams as indicated by the high N(P) for 
P3 > 0.9 (blue histogram in Figure [5J . Caution is there- 
fore advisable in the evaluation of individual ultra-faint 
systems given their propensity to appear physically close 
to the Sgr streams simply by virtue of the location of the 
SDSS survey fields. 

4.2. Caveats 

Our numerical method of quantifying the degree of as- 
sociation of individual Milky Way satellites with the Sgr 
stellar debris streams is based upon some significant as- 
sumptions. First, we assume that the LM10 model is an 
accurate representation of the real Sgr stream. While 
this is a good assumption for the L1Y and T1Y wraps 
of the leading and trailing tidal arms (where the model 
was constrained to match observational data), it is un- 
known how well the model may describe tidal debris more 
widely separated from the dwarf. It may be the case that 
the most extended regions of the tidal tails (if they exist 
in the Galactic halo) depart significantly from the LM10 
model. Indeed, since the LM10 model does not explicitly 
include orbital evolution of the Sgr dwarf (e.g., via dy- 
namical friction), some mismatch in the oldest sections 
of the tidal streams is likely. We therefore generally place 
our greatest confidence in the P3 statistic, which quanti- 
fies the association of individual systems with only those 
segments of the Sgr stream that have been conclusively 
detected in the Galactic halo and used to constrain the 

8 We use an SDSS DR7 lookup table kindly provided by M. Juric 
2010 (priv. comm.). 



LMIO model. 

Second, we assume that globular clusters stripped from 
Sgr follow the path of Sgr tidal debris. We believe this 
to be a reasonable assumption because neither unbound 
tidal debris nor globular clusters are sufficiently massive 
to have their orbits noticeably altered by dynamical fric- 
tion in the gravitational potential of the Milky Way. 

Third, and particularly relevant to the case of the 
ultra-faint satellites, we have assumed that the artificial 
comparison satellite distribution is spherically symmetric 
about the Milky Way. This assumption, and the ensu- 
ing P statistics derived from it, may be incorrect if the 
real satellite population is distributed in a strongly non- 
isotropic manner. Such a strongly flattened distribution 
is not observed for the Galactic globular cluster popu- 
lation, but may prove to be the case for the ultra-faint 
satellite population if they tend to be aligned with the 
observed plane of Milky Way satellite galaxies (see, e.g., 
Zentner et al. 2005; Metz et al. 2009), or if they fell into 
the Galactic halo in a group with the Sgr dwarf (see, e.g., 
Li & Helmi 2008; D'Onghia & Lake 2009). 

Finally, we assume that no Sgr clusters have been 
tidally destroyed over the past few Gyr. While this would 
not affect our identification of which current satellites 
originated in the Sgr system, it may cause us to under- 
estimate the total globular cluster content of the origi- 
nal system. For clusters still contained within the Sgr 
dwarf, this is likely a reasonable assumption (see, e.g., 
Penarrubia et al. 2009). Once unbound from the Sgr 
gravitational potential however, the validity of this as- 
sumption is unknown since some globular clusters are 
known to be tidally disrupting in the Galactic potential 
(e.g., Pal 5; Odenkirchen et al. 2001), even at significant 
distances from the Galactic center. 

5. ANALYSIS OF INDIVIDUAL STAR CLUSTERS 

In Figures |5] and [7] we overplot each of the satellites 
from Table [2] (coded by ID number) against the LM10 
model of angular positions, distances, and radial veloci- 
ties for the Sgr trailing/leading arm streams respectively. 
For comparison, in Figure |8] we overplot the satellites 
against 2MASS M-giants whose radial velocities (Law et 
al. 2004; Majewski et al. 2004) are consistent with mem- 
bership in the Sgr stream (see Figure 1 of LM10). Satel- 
lites that we consider are/are-not candidates for mem- 
bership in the Sgr stream based on the selection criteria 
P3 < 0.15 or Pj < 0.15 are summarized in Tables [3] and 
0] respectivelyO Since P3 best corresponds to the obser- 
vationally verified sections of the Sgr stream, we focus 
our discussion primarily on this statistic, except in those 
cases when a better match is found using Pj. 

Given the large number of globular clusters (7 ± 2) 
that we expect to fall in the range P < 0.15 by virtue of 
chance superposition in angular position, distance, and 
radial velocity with the Sgr stream, not all candidates 
listed in Table [3] will prove to be genuinely associated 
with Sgr. Instead P is a useful statistical indicator, sug- 
gesting that a system is sufficiently well-matched to war- 
rant consideration; these objects are discussed further in 

9 The four globular clusters AM 4, IC 4499, BH 176, and ESO 
280-SC06, and the ultrafaint dwarf Segue 3 do not presently have 
published radial velocities. We therefore estimate P3 and P-j for 
these objects using a revised version of the \ statistic which is 
based upon only distance and angular position. 
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ID 


Name 


Pa 


Pr 


^~3,T1Y 


^3,L1Y 


•FV.Tl 


Pr,T2 


Pr,u 


^7,L2 










Star 


Clusters 










1 


NGC 288 


0.095 


0.315 


0.095 


0.934 


0.315 


0.727 


0.963 


0.366 


3 


Whiting 1 


0.001 


0.010 


0.001 


0.796 


0.010 


0.995 


0.588 


0.064 


7 


Pal 2 


0.123 


0.316 


0.123 


0.971 


0.316 


0.946 


0.327 


0.659 


11 


Berkeley 29 


0.027 


0.039 


0.378 


0.027 


0.750 


0.743 


0.039 


0.608 


IS 


NGC 4147 


0.058 


0.132 


0.058 


0.317 


0.132 


0.999 


0.702 


0.240 


21 


M 53 


0.106 


0.112 


0.106 


0.185 


0.147 


0.805 


0.422 


0.112 


22 


NGC 5053 


0.040 


0.032 


0.040 


0.274 


0.032 


0.999 


0.591 


0.442 


23 


M 3 


0.344 


0.122 


0.412 


0.344 


0.412 


0.732 


0.690 


0.122 


27 


NGC 5634 


0.004 


0.031 


0.004 


0.535 


0.031 


0.589 


0.712 


0.608 


31 


Pal 5 


0.045 


0.117 


0.045 


0.202 


0.117 


0.582 


0.499 


0.853 


.37 


Pal 15 


0.128 


0.390 


0.890 


0.128 


0.982 


0.823 


0.390 


0.999 


39 


IC 1257 


0.695 


0.112 


0.695 


0.698 


0.761 


0.990 


0.962 


0.112 


42 


M 54 


0.001 


0.001 


0.001 


0.975 


0.001 


0.794 


0.410 


0.514 


-11 


Terzan 7 


0.122 


0.090 


0.122 


0.869 


0.090 


0.764 


0.549 


0.458 


45 


Arp 2 


0.001 


0.015 


0.001 


0.615 


0.015 


0.660 


0.344 


0.277 


16 


Terzan 8 


0.001 


0.009 


0.001 


0.479 


0.009 


0.654 


0.441 


0.640 


53 


Pal 12 


0.012 


0.040 


0.012 


0.048 


0.066 


0.865 


0.040 


0.835 


Ultrafaint D warts 


58 


Segue 1 


0.318 


0.113 


0.318 


0.999 


0.113 


0.970 


0.999 


0.483 


61 


Bootes II 


0.124 


0.239 


0.918 


0.124 


0.566 


0.662 


0.239 


0.758 



a Probability P that a randomly selected artificial satellite would match the indicated segment of the Sgr stream as well as, or better than, the 
indicated satellite. Values P < 0.15 arc indicated in bold font. 



the subsections below. In many cases, analysis of proper 
motions or other physical considerations can help confirm 
whether an association is likely to be genuine. Where 
available, proper motions for each of the satellites listed 
in Table [5] have been overplotted against those of the 
model Sgr stream in Figures 1^1 and ITU1 for trailing/leading 
arms respectively. In the interests of brevity, we do not 
discuss satellites with both P3 > 0.15 and P-j > 0.15 
(which are extremely unlikely to belong to Sgr) in detail 
unless they exhibit some other point of interest, or have 
been described as Sgr cluster candidates previously in 
the literature. 

In the following subsections we begin our discussion 
with the classical Sgr 'core' globular clusters M 54, Arp 
2, Terzan 7, and Terzan 8, and thereafter proceed to dis- 
cuss individual systems in roughly alphanumerical order 
to aid the reader in easily locating particular sections of 
interest. In most cases, we give the A Q coordinate and 
ID number of each cluster in the corresponding subsec- 
tion heading. We group our discussion of the clusters 
M2, NGC 5 466, NGC 5824, NGC 6426, and Rup 106 
together in £15.171 and discuss the four clusters without 
radial ve locity data (AM 4, BH 176, ESO 280-SC06, IC 
4499) in ^5.181 Finally, we also discuss globular clusters 
that did not meet our initial 8 < rcc < 100 kpc selection 
criterion briefly in §5.191 Our conclusions regarding the 
star cluster budget of Sgr are summarized in Section [7] 

5.1. M 54, Arp 2, Terzan 7 and Terzan 8 

The four globular clusters M 54, Arp 2, Terzan 7, and 
Terzan 8 have traditionally been associated with the Sgr 
dSph (e.g., Ibata et al. 1994; Da Costa & Armandroff 
1995); these four 'core' clusters lie within ~ 1 kpc of 
the best-fitting Sgr debris plane and 10° (~ 5 kpc) of 
the Sgr dwarf core along the trailing tidal stream. To 
within observational uncertainty, M 54 appears to be 
spatially coincident with the core of Sgr, and to lie at 
the same distance and radial velocity. While M 54 is cer- 
tainly associated with the Sgr dwarf (we calculate a less 



than P3 = 0.1% chance that an artificial globular cluster 
would chance to match Sgr so well), its velocity disper- 
sion and luminosity profile indicate that it formed inde- 
pendently of the Sgr core and later sank to its present 
position via dynamical friction (Bellazzini et al. 2008, 
see also Majewski et al. 2003; Monaco et al. 2005; Siegel 
et al. 2007). 

Arp 2 and Terzan 8 each have similarly strong probabil- 
ities for association. While proper motion information is 
not available for these clusters to confirm the association, 
the strength of the match in angular position, distance, 
and radial velocity is sufficient to give us confidence in 
concurring with previous studies (see Table Q} that M 54, 
Arp 2, and Terzan 8 are genuine Sgr members. Indeed, 
out of a sample of 51 randomly chosen artificial globular 
clusters we should expect less than 0.05 of them to match 
the Sgr stream so well by chance, let alone the 3 that are 
actually observed. 

In contrast, Terzan 7 has a significantly lower confi- 
dence for association (P3 rs 0.12); we expect ~ 6 satel- 
lites in a given sample of 51 to meet this threshold sim- 
ply by chance. However, the high P3 primarily reflects 
our adopted distance to Terzan 7 (d = 23.2 kpc), which 
is almost 5 kpc less than the distance to the Sgr core 
(£>sgr = 28 kpc; Siegel et al. 2007) adopted in the 
LM10 model. If the LM10 model had instead adopted 
a slightly smaller distance for Sgr (reasonable estimates 
range from Dg gI — 24 — 28 kpc; see Table 2 of Kun- 
der & Chaboyer 2009 for a summary), the quality of fit 
of Terzan 7 to the Sgr stream would improve consider- 
ably. Similarly, had we adopted a distance of d = 25.75 
kpc for Terzan 7 (M. Siegel 2009, private communica- 
tion) based on new photometry akin to that used to set 
the Sgr distance scale in Siegel et al. (2007), we would 
have found P 3 = 0.032 and P 7 = 0.052. While Terzan 7 
is therefore not as robust a Sgr cluster candidate as M 
54, Arp 2, or Terzan 8, physical association is still likely. 
Future proper motion measurements should be able to 
verify this; the proper motion of the Terzan 7 region of 
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TABLE 4 
Systems Unlikely to Match the Sgr Stream 51 



ID 


Name 


P 3 


Pi 


P?,,T1Y 


^"3,L1Y 


^V,T1 


Pt,T2 


*V,Li 


^7,L2 


Star Clusters 


2 


NGC 362 


0.463 


0.403 


0.463 


0.715 


0.841 


0.952 


0.403 


0.910 


1 


NGC 1261 


0.233 


0.371 


0.233 


0.948 


0.602 


0.767 


0.982 


0.371 


5 


Pal 1 


0.773 


0.754 


0.842 


0.773 


0.962 


0.971 


0.912 


0.754 


6 


Eridanus 


0.495 


0.482 


0.495 


0.999 


0.500 


0.999 


0.755 


0.482 


8 


NGC 1851 


0.708 


0.552 


0.834 


0.708 


0.996 


0.820 


0.552 


0.736 


9 


M 79 


0.356 


0.395 


0.356 


0.416 


0.755 


0.901 


0.395 


0.498 


10 


NGC 2298 


0.642 


0.803 


0.809 


0.642 


0.990 


0.999 


0.803 


0.900 


12 


Saurer 1 


0.190 


0.313 


0.677 


0.190 


0.851 


0.773 


0.313 


0.385 


13 


NGC 2419 


0.386 


0.487 


0.386 


0.999 


0.671 


0.487 


0.673 


0.729 


11 


Pyxis 


0.612 


0.722 


0.999 


0.612 


0.986 


0.856 


0.722 


0.769 


15 


NGC 2808 


0.962 


0.992 


0.999 


0.962 


0.999 


0.999 


0.999 


0.992 


16 


Pal 3 


0.825 


0.540 


0.953 


0.825 


0.842 


0.646 


0.540 


0.639 


17 


NGC 3201 


0.800 


0.937 


0.882 


0.800 


0.937 


0.999 


0.986 


0.947 


19 


Rup 106 


0.615 


0.781 


0.615 


0.763 


0.781 


0.938 


0.982 


0.782 


20 


M 68 


0.608 


0.266 


0.608 


0.739 


0.744 


0.856 


0.970 


0.266 


24 


NGC 5286 


0.757 


0.878 


0.757 


0.856 


0.878 


0.964 


0.993 


0.954 


25 


AM 4 b 


0.367 


0.597 


0.437 


0.367 


0.605 


0.774 


0.597 


0.612 


26 


NGC 5466 


0.531 


0.505 


0.531 


0.538 


0.505 


0.999 


0.853 


0.938 


28 


NGC 5694 


0.414 


0.385 


0.414 


0.717 


0.387 


0.639 


0.970 


0.385 


29 


IC 4499 b 


0.980 


0.885 


0.980 


0.999 


0.999 


0.885 


0.992 


0.999 


30 


NGC 5824 


0.287 


0.438 


0.287 


0.579 


0.438 


0.731 


0.876 


0.837 


32 


BH 176 b 


0.828 


0.845 


0.933 


0.828 


0.989 


0.999 


0.979 


0.845 


33 


Pal 14 


0.353 


0.572 


0.968 


0.353 


0.935 


0.919 


0.572 


0.999 


34 


NGC 6101 


0.978 


0.931 


0.999 


0.978 


0.967 


0.931 


0.999 


0.961 


35 


M 13 


0.628 


0.689 


0.729 


0.628 


0.689 


0.798 


0.939 


0.849 


36 


NGC 6229 


0.424 


0.502 


0.424 


0.490 


0.502 


0.956 


0.777 


0.926 


38 


M 92 


0.551 


0.901 


0.782 


0.551 


0.901 


0.972 


0.901 


0.991 


40 


NGC 6426 


0.804 


0.364 


0.877 


0.804 


0.952 


0.999 


0.988 


0.364 


-11 


ESO 280-SC06 b 


0.872 


0.453 


0.882 


0.872 


0.453 


0.944 


0.982 


0.760 


43 


M 56 


0.999 


0.966 


0.999 


0.999 


0.972 


0.999 


0.999 


0.966 


47 


M 75 


0.346 


0.474 


0.862 


0.346 


0.999 


0.523 


0.648 


0.474 


48 


NGC 6934 


0.795 


0.622 


0.915 


0.795 


0.999 


0.897 


0.622 


0.990 


49 


M 72 


0.647 


0.444 


0.875 


0.647 


0.999 


0.854 


0.444 


0.945 


50 


NGC 7006 


0.841 


0.535 


0.854 


0.841 


0.996 


0.761 


0.535 


0.944 


51 


M 15 


0.425 


0.278 


0.425 


0.548 


0.806 


0.956 


0.278 


0.903 


52 


M 2 


0.333 


0.181 


0.338 


0.333 


0.719 


0.955 


0.181 


0.986 


51 


Pal 13 


0.366 


0.155 


0.772 


0.366 


0.986 


0.711 


0.155 


0.982 


55 


NGC 7492 


0.318 


0.717 


0.318 


0.929 


0.717 


0.760 


0.962 


0.788 










Ultrafaint Dwai 


■Is 










56 


Segue 2 


0.892 


0.330 


0.892 


0.999 


0.999 


0.948 


0.716 


0.330 


57 


UMall 


0.954 


0.185 


0.971 


0.954 


0.999 


0.185 


0.956 


0.927 


59 


Willman I 


0.940 


0.694 


0.940 


0.970 


0.694 


0.707 


0.999 


0.887 


60 


Coma Ber 


0.446 


0.434 


0.791 


0.446 


0.434 


0.999 


0.791 


0.818 


62 


Bootes I 


0.493 


0.774 


0.999 


0.493 


0.878 


0.999 


0.774 


0.983 


63 


Segue 3 b 


0.945 


0.833 


0.945 


0.968 


0.985 


0.999 


0.833 


0.938 


64 


Pisces I 


0.999 


0.290 


0.999 


0.999 


0.950 


0.290 


0.986 


0.935 



the Sgr stream is extremely well constrained by current 
observational data, and if Terzan 7 is moving with the 
stream the LM10 model suggests that we should expect 
it to have a proper motion fi a cosd = —2.4 ± 0.1 mas 
yr _1 and [i$ = —1.4 ± 0.1 mas yr _1 . 



a Probability P that a randomly selected artificial satellite would match the indicated segment of the Sgr stream as well as, or better than, the 
indicated satellite. Values P < 0.15 are indicated in bold font. 

No radial velocity information was available, so association statistics are calculated from angular position and distance alone. 

the Sgr stream{3 Although typically regarded as open 
clusters, Be 29 and Sa 1 have ages, metallicities, and 
richness intermediate between classical globular cluster 
and open cluster designations. Since Be 29 in particular 
(1) lies ~ 2 kpc above the Galactic plane, (2) is the most 
distant known open cluster (Tosi et al. 2004), and (3) is 
brighter and richer than some globular clusters (e.g., AM 
4, Whiting 1, and Fornax-Hl), we consider it reasonable 
that these two clusters may indeed be part of the halo 
population rather than the disk. 

While both clusters are approximately coincident with 
the Sgr LI arm, Be 29 is the superior match with P3 = 
0.027 (compared to P 3 = 0.190 for Sa 1). Visually, Figure 

10 Note that the declination given for Be 29, and the right ascen- 
sion and declination given for Sa 1 in Table 1 of Carraro & Bensby 



5.2. Berkeley 29 and Saurer 1 (ID=11/12, 

A Q = 177 /182 J 

Previous studies of the star cluster budget of the Sgr 
dwarf (e.g., Palma et al. 2002; Bellazzini et al. 2003a) 
have typically neglected sparse open clusters located in 
the plane of the Galactic disk. However, Carrarro & 
Bensby (2009) argue that the metal-poor ([Fe/H] = 
-0.44 ± 0.18 and -0.38 ± 0.14 respectively; Carraro et 
al. 2004) cluster pair Be 29 and Sa 1 may be members of 



(2009) are incorrect. 
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Fig. 6. — Plot shows the predicted angular position Bq, distance d, and radial velocity i>gsr as a function of Sgr orbital longitude Am 
for the LM10 model of the Sgr trailing tidal stream. Different colors for points represent debris lost on different orbital passages (see Sj2]l. 
Milky Way stellar subsystems are overplotted, ID numbers correspond to those defined in Table [2] Typically the stellar system is located 
at the bottom left corner of the corresponding ID label, although the locations of some ID numbers have been adjusted slightly to minimize 
confusion from overlapping labels. 



[7] confirms this: Be 29 is superimposed on the LI arm 
while Sa 1 lies off the projected path of the Sgr stream 
by ~ 20° (4.1cr) in Sgr orbital latitude Bq. Both clusters 
lie at a distance d = 13.2 kpc consistent with that of 
the LI arm, but Be 29 is a slightly better match to the 
LI stream in radial velocity (although still discrepant by 
~ 60 km s _1 , or 3<r) than is Sa 1 (off by ~ 80 km s _1 , 
or 4cr). The radial velocity of the LI stream is not well 
constrained in this region however since it lies ~ 50° and 
~ 20° respectively beyond the last 2MASS (see, e.g., 
Figure H]) and SDSS observations used to constrain the 
LM10 model, and it is therefore possible that the LM10 
model at the location of Be 29 and Sa 1 may be in error 
by up to a few tens of kilometers per second. 

The simplest hypothesis is that neither Sa 1 nor Be 29 
are associated with the Sgr stream, and that both are 
unusual Galactic disk open clusters that happen to lie 
near the nexus at which the LI arm crosses through the 
disk and happen to have radial velocities resembling that 



of the stream. Given the thousands of open clusters in 
the Galactic disk, it is perhaps unsurprising that some 
cluster might happen to have the requisite combination 
of properties. This simple hypothesis is probably true 
for Sa 1 considering its low statistical probability for as- 
sociation with the Sgr stream (P3 = 0.190)E3but the low 
P3 = 0.027 for Be 29 may favor the Sgr-origin hypothesis. 
Although a tentative measurement of the proper motion 
of Be 29 was presented by Dias et al. (2006), Frinchaboy 
(2006) demonstrated that the majority of stars used to 
make this measurement were foreground field stars. Only 
a single star had the correct magnitude and radial veloc- 
ity to be a cluster member, and the proper motion of this 
star (/j, a cosS = — 2.0±7.5 mas yr -1 , (is = — 5.7±7.5 mas 
yr _1 ; P. Frinchaboy 2010 [priv. communication]) is too 

11 Alternatively it may instead be associated with some other 
tidal feature such as Monoceros/GASS; see Frinchaboy et al. 2006, 
Pcfiarrubia et al. 2005. 
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Fig. 7. — As FigureO but for the leading tidal stream. 



uncertain to offer a significant constraint on the relation 
between Be 29 and the Sgr stream. 

If Be 29 is indeed a member of the Sgr stream, it is in- 
triguing that this cluster is wrapped farther from the Sgr 
dwarf than all tidal debris conclusively observed to date, 
and is the only known cluster (except possibly for Pal 12) 
associated with the leading tidal arm. The LI arm at the 
location of Be 29 has been unbound from the Sgr core 
for at least 2.5 Gyr, meaning that if Be 29 formed in the 
Sgr dwarf ~ 4 Gyr ago (Carraro et al. 2004 derive an age 
of 4.5 Gyr, while Tosi et al. 2004 find ~ 3.4 - 3.7 Gyr) 
it must have been torn from Sgr shortly after its forma- 
tion, and has since experienced at least 3 orbits around 
the Milky Way. Given the similarities between Be 29, 
Whiting 1, and Terzan 7 in age, metallicity, luminosity, 
horizontal-branch type, and size (see discussion in Sj7|) we 
suggest that Be 29 may be a paradigm for the origins of 
some old, metal-poor 'open clusters' in the Milky Way; 
namely, that they may be a class of sparse globular clus- 
ter formed in dSph, stripped from their sheltering birth 
galaxies, and subjected to the gravitational tides of the 
Milky Way. 



5.3. IC 1257 (ID=39, A Q = 334°; 

IC 1257 has never previously been suggested as a pos- 
sible Sgr cluster, but our analysis indicates that it re- 
sembles the L2 wrap at A = 334° to P 7 = 0.112. This 
association is not particularly compelling however; while 
it closely matches the distance of the L2 wrap (~ 25 kpc) 
the angular position is inconsistent by ~ 10° and the ra- 
dial velocity by ~ 115 km s _1 . It is only due to the 
large apparent width of the L2 stream in radial velocity 
(o"v,Sgr ~ 60 km s _1 ) that the association statistic is cal- 
culated to be as good as P7 = 0.112. Since P7 = 0.112 
corresponds to ~ 6 false associations out of a sample of 
51 clusters, and there is not even any concrete evidence 
that the L2 wrap actually exists in the Galactic halo, we 
doubt that IC 1257 is genuinely associated with the Sgr 
stream. 

5.4. M 3 (NGC 5272; ID=23, A Q = 265°; 

At a distance of 10.4 kpc, M 3 is too close to be as- 
sociated with anything except the secondary L2 wrap of 
the leading arm. As illustrated in Figures [7] and [8j the 
radial velocity of M 3 is consistent with the L2 identifi- 
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Fig. 8. — As Figures \E\ and but black solid points indicate observed Sgr stream M-giants from the leading and trailing arm samples 
described by LM10. 



cation, but the angular position is discrepant by 26.5°. 
The importance of this angular discrepancy is mitigated 
somewhat by the large angular width of the L2 stream 
(26.5° corresponds to 3.8a), but there is a Pj = 12.2% 
chance (i.e., 6 out of a sample of 51) that a randomly cho- 
sen artificial cluster would appear to be associated with 
some wrap of the Sgr stream more convincingly than is 
M 3. 

Adopting the proper motions indicated in Table [21 we 
observe (Figure IT0|) that although the us for M 3 is con- 
sistent with the L2 wrap of the Sgr stream, /i Q cos<5 is 
inconsistent by ~ 6.9 mas yr _1 (corresponding to 2.3c). 
While at least four different values have been measured 
for the proper motion of M 3 (see compilation in Palma et 
al. 2002), none of these measurements agree more closely 
with the proper motion expected for the Sgr stream, and 
some (e.g., Cudworth & Hansen 1993) are significantly 
poorer matches to the estimated /j,g. We therefore con- 
clude that M 3 is probably not physically associated with 
the Sgr stream, although future constraints on the Sgr 
stream at extremely large separations from the dSph, and 
measurements of the proper motion of M 3 to an accu- 



racy ~ 0.5 mas yr 1 will help to check this conclusion. 

5.5. M 53 (NGC 5024; ID =21, A© = 265° J 

M 53 is situated towards the North Galactic Cap and 
has been previously discussed as a possible candidate Sgr 
globular cluster by Palma et al. (2002; although these 
authors rule it out on the basis of its low metal content 
[Fe/H] = -1.99), Bellazzini et al. (2003a), and Martinez 
Delgado et al. (2004). We find a mild probability for 
association between M 53 and the Tl (P3 = 0.106), LI 
(P3 = 0.185; barely failing the P < 0.15 criterion), and 
L2 {P-j = 0.112) arms. Figures [6] and [7] indicate that 
while M 53 matches the Tl arm in angular position and 
distance it differs by 96 km s _1 (i.e., 4c) in radial velocity, 
and while it matches the L1/L2 arms in radial velocity 
it lies at the edge of both streams and is simultaneously 
too close for the LI arm and too distant for the L2 arm 
(17.8 kpc versus 45 ± 4 kpc and 11 ± 3 kpc respectively). 

While the proper motion (Figures [5] and [TU)) along 
the declination axis (/is) is consistent with all three 
arms to within observational uncertainty, the proper 
motion along right ascension (/i Q cosd>) is inconsistent 
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Fig. 9. — Plots predicted proper motion along right ascension (/i Q cos<5) and declination (fJ-s) a s a function of Sgr orbital longitude Aq 
for the LM10 model of the Sgr trailing tidal stream. Different colors for points represent debris lost on different orbital passages (see 
$2^ . Milky Way stellar subsystems are overplotted, ID numbers correspond to those defined in Table [2] Errorbars on the proper motion 
measurements represent la. The locations of ID numbers have been adjusted to minimize confusion from overlapping labels. 



by 4.7/1.7/10.2 mas yr" 1 (A.3a/1.7a/i.la) with the 
T1/L1/L2 arms respectively. Given these difficulties in 
associating M 53 conclusively with any particular arm of 
the Sgr stream and the likelihood that ~ 5 clusters with 
similarly good P statistics will be found in our sample of 
51 clusters by chance, we conclude that M 53 is probably 
not physically associated with the Sgr system. 

5.6. NGC 288 (ID=1, A Q = 76°; 

NGC 288 lies within ~ 1 kpc of the Sgr plane in the 
direction of the South Galactic Pole, and (as pointed out 
by Bellazzini et al. 2003a) has a radial velocity very close 
to that expected for the Tl trailing arm. Although the 
distance to NGC 288 (d = 8.8 kpc) is consistent with 
that of the M-giants (Figure [SJ), this is likely due to the 
large observational uncertainty in the photometric paral- 
laxes of the M-giants; NGC 288 lies ~ 10 kpc closer than 
the LM10 model of the Tl stream. Since the association 
probability P3 = 0.095 corresponds to ~ 5 false positives 
anticipated in our sample of 51 globular clusters, we are 



inclined to doubt the association of NGC 288 with the 
Sgr stream. Table [3] indicates that significant matches 
to other segments of the Sgr stream are even less likely 
This negative hypothesis is borne out by the proper mo- 
tions: Figure [9] demonstrates that the \x a cose) and /is 
of NGC 288 are discrepant with the Tl stream at this 
location by 5.9 mas yr _1 and 2.4 mas yr _1 respectively. 
Combining the uncertainty in the observed proper mo- 
tions and the variance in proper motions of particles in 
the LM10 model, this corresponds to an offset of ~ 22c. 
We are therefore confident in our conclusion that, despite 
its proximity to the Sgr plane, NGC 288 is not physically 
associated with the Sgr stream. 



5.7. 



NGC 2419 (ID=13, A = 190°; 

Irwin (1999) and Newberg et al. (2003) note that the 
distant (d = 84 kpc) Galactic globular cluster NGC 2419 
([l,b] = [180°, 25°]) lies near the Sgr debris plane, and 
suggest that it may be associated with the Sgr trailing 
tidal arm. Indeed, Figure [6] shows that NGC 2419 lies 
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Fig. 10. — As Figure [9] but for the leading tidal stream. 



in projection against the Tl arm with B & « —9°. How- 
ever, the large distance to the cluster is poorly matched 
by the LM10 model, which predicts that the Tl stream 
reaches apogalacticon at a distance ~ 61 kpc around at 
A© ~ 170°, and approaches closer to the Milky Way 
again at greater orbital longitudes. NGC 2419 is for- 
mally inconsistent with the LM10 model by ~ 38 kpc in 
distance and ~ 160 km s _1 in radial velocity, giving a 
poor association statistic P3 = 0.386. 

Direct examination of the Tl M-giant stream observed 
by Majewski et al. (2003; their Figure 10) is less clear; 
while it admits the possibility that the observed stream 
may have a distance ~ 80 kpc at A = 190° (perhaps 
tracing a previously larger orbit for Sgr, which has since 
sunk deeper into the Galactic gravitational potential via 
dynamical friction), the stream at this location is diffi- 
culty to trace because it is poorly populated by M-giants 
and has an uncertain distance calibration. While the 
LM10 model of the Tl arm is not directly constrained at 
A Q = 190°, it agrees well with observations of this arm 
at both A Q - 160° (i.e., the M-giants) and A - 225° 
(see discussion in H6.5I) , suggesting that it should not be 



wildly inaccurate in the longitude range of NGC 2419. 
We therefore conclude that NGC 2419 is unlikely to be 
genuinely associated with the Sgr dSph. 

5.8. NGC 4147 (ID=18, A = 251°; 

NGC 4147 (1.4 kpc from the Sgr midplane) has been 
previously discussed as a potential member of the Sgr 
system by Bellazzini et al. (2003ab) and Palma et 
al. (2002). This cluster is a formally fair match to 
the Tl principle wrap of the trailing tidal arm (as ob- 
served in the direction of the North Galactic Cap) , with 
P 3 = 0.058. As indicated by Figure[6]however, NGC 4147 
lies near the edge of the Tl stream in projection and is 
discrepant from its radial velocity by ~ 86 km s _1 . NGC 
4147 is therefore not an extremely attractive candidate 
for membership in the Sgr stream because ~ 3 such clus- 
ters in our sample of 51 are expected to appear at least 
as closely associated by chance. 

While proper motions have been measured for NGC 
4147 (Brosche et al. 1991), they do not conclusively con- 
firm or deny membership in the Sgr stream; as illustrated 
in Figure [9] \x a cos<5 is discrepant from the Tl stream by 
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2.3 mas yr _1 (corresponding to ~ 1.5a), while fi$ is con- 
sistent with the Tl stream. However, the uncertainties 
on the observed proper motions (particularly for [i$) are 
so large that they contain little discriminatory power. 
We are therefore unable to determine conclusively the 
whether NGC 4147 is a member of the Sgr stream, al- 
though it appears to be only a weak candidate. 

5.9. NGC 5053 (ID =22, A Q = 266°; 

NGC 5053 lies near M 53 in the North Galactic Cap 
(separation ~ 1°) and at a similar distance (as discussed 
by Palma et al. 2002; Bellazzini et al. 2003a), although 
the two clusters have radial velocities that differ by ~ 125 
km s _1 . As illustrated by Figure^ the radial velocity of 
NGC 5053 (ugsr = 34.1 ± 0.4 km s" 1 ) is a much better 
match to that expected for the Tl stream (fGSR = 25±17 
km s ) than was M 53. This results in a much improved 
value P3 = 0.04; only ~ 2 such systems are expected to 
occur by chance in our sample of globular clusters. 

Lauchner et al. (2006; see also Chun et al. 2010) have 
reported a tentative detection of tidal streams emanating 
from the cluster, stretching from (a, 8) = (199°, 18°) to- 
wards (195°, 15°). If genuine, these streams may weaken 
the case for association of NGC 5053 with the Sgr stream 
since they suggest a motion oriented at a transverse an- 
gle ~ 60° to the Sgr plane. However, these streams may 
be too close to the cluster (^3 — 4 tidal radii; Lauchner 
ct al. 2006) to provide a good indication of the actual 
direction of the cluster's orbit (see, e.g., Johnston et al. 
2002; Montuori et al. 2007; Klimentowski et al. 2009; 
Odenkirchen et al. 2009), and cannot yet be used to 
eliminate NGC 5053 from the Sgr family. 

Unfortunately, no proper motion has been measured 
for NGC 5053, and it is therefore not possible to conclu- 
sively determine its membership status in the Sgr stream. 
If it proves to be a genuine member, we might expect the 
proper motion to be similar to \i a cos(5 = —4.2 ± 0.5 mas 
yr" 1 and fi$ = —0.1 ±0.1 mas yr _1 . At present however, 
the case for membership appears to be fairly good. If 
confirmed, NGC 5053 would be the most metal poor of 
the Sgr globular clusters with [Fc/H] = —1.98 (Marin- 
Franch et al. 2009). 

5.10. NGC 5634 (ID=27, A = 293° ) 

NGC 5634 lies within 0.5 kpc of the Sgr plane, and is 
consistent with the Tl segment of the trailing arm stream 
which has been wrapped around the Milky Way into the 
direction of the North Galactic Cap. As indicated by 
Tabled there is a probability P3 = 0.4% for a randomly 
chosen halo globular cluster to match the Sgr stream this 
well, corresponding to ~ 0.2 such clusters in a given sam- 
ple of 51. As illustrated in Figure [51 NGC 5634 matches 
the model Tl stream to within ~ 6° (la) in B & , 2 kpc 
(< la) in distance, and 4 km s _1 (< la) in i>gsr- The 
extremely low metallicity of NGC 5634 ([Fe/H] = -1.94; 
Bellazzini et al. 2002) is consistent with that of the three 
Sgr core clusters M 54, Arp 2, and Terzan 8, which have 
-1.2 > [Fe/H] > -1.8. We agree with the possibility 
suggested by Bellazzini et al. (2002; 2003a) that NGC 
5634 is an extremely strong candidate for a Sgr globular 
cluster. 

Proper motion measurements have not yet been ob- 
tained to test this hypothesis, but if NGC 5634 is indeed 



associated with the Tl stream we might expect it to have 
Ha cos<5 k, —3.8 ± 0.4 mas yr _1 and /j,s « +0.2 ±0.1 mas 
yr _1 . 

5.11. Pal 2 (ID=7, A = 154°; 

Pal 2 is a heavily obscured (E(B - V) = 1.24 ± 0.07; 
Harris et al. 1997) globular cluster that lies in the di- 
rection of the Galactic anticenter. The possibility that 
Pal 2 may be part of the Tl trailing arm of Sgr has been 
previously noted by Lynden-Bell & Lynden-Bell (1995), 
Irwin (1999), and Majewski et al. (2004). 

As illustrated in Figure [51 Pal 2 lies in the Tl trail- 
ing arm stream in projection, and has a radial veloc- 
ity consistent with this identification, but at a distance 
d = 28 ± 4 kpc (Harris et al. 1997), roughly half that ex- 
pected for the Tl stream at this orbital longitude in the 
LM10 model (d = 55 ±4 kpc). This 27 kpc distance mis- 
match translates to a P3 m 12% chance that a random 
halo globular cluster would match the Sgr stream as well 
as, or better than, Pal 2 by chance. Since we expect six 
such false positives in our sample of 51 clusters, we are 
inclined to treat the potential association of Pal 2 with 
the Sgr stream with skepticism. 

It is possible since Pal 2 is so heavily obscured that its 
distance (Harris et al. 1997) may have been underesti- 
mated, although we note that Sarajedini et al. (2007) 
derive a similar distance d = 27 ± 2 kpc based on main- 
sequence fitting. Alternatively, the LM10 model may sys- 
tematically overestimate the distance to the Sgr stream 
at this location. As discussed by Majewski et al. (2004), 
Pal 2 appears to lie only a few kpc closer than the clos- 
est of the M giants observed in this region (Majewski 
et al. 2003). However, the M-giants distances in this 
region may be unreliable given the evolving metallicity 
distribution function (MDF) along the Sgr tidal streams 
(Chou et al. 2007, 2010); the LM10 model was designed 
to account for the evolution in the MDF and matches the 
observed M-giants in apparent magnitude (see Figure 18 
of LM10). Unfortunately, the proper motion of Pal 2 has 
not been measured and it is therefore not presently pos- 
sible to determine conclusively the relation of Pal 2 to 
the Sgr stream. 

5.12. Pal 5 (ID '=31, A© = 300°; 

Pal 5 lies near NGC 5634, and within ~ 4 kpc of the 
Sgr plane. At a distance of 23.2 kpc it falls inside the 
loop of Sgr LI leading arm debris, and is consistent with 
the distance expected for the Tl stream wrapped into 
the northern Galactic hemisphere (see Figure [5J . How- 
ever, the radial velocity of Pal 5 is inconsistent with this 
identification by ~ 83 km s _1 (3.8er). This radial veloc- 
ity mismatch gives a probability statistic P3 = 0.045, 
corresponding to ~ 2 such false positives expected in our 
sample of 51 globular clusters. 

The situation of Pal 5's possible association with Sgr 
has been discussed by Bellazzini et al. (2003a) based 
on its close proximity to the Sgr plane, by Palma et al. 
(2002) on the basis of the location of its orbital pole, and 
with regard to its proper motion by Scholz et al. (1998). 
It has also been noted (Palma et al. 2002) that Pal 5 
has a similar mass and concentration to the candidate 
Sgr globular clusters Arp 2, Terzan 7, Terzan 8, and Pal 
12, and that its alpha element abundance is similar to 
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that observed in M 54, Ter 7, and Pal 12 (Sbordone et 
al. 2005). While these latter authors suggest that Pal 5 
may be a member of the Sgr stream based on its chemical 
properties, they note that its orbital characteristics sug- 
gest it is more likely a candidate to be associated with 
older rather than younger tidal debris phases. Indeed, 
Palma et al. (2002) demonstrate that while the orbital 
poles of Sgr and Pal 5 are somewhat similar (separated 
by about 27°), the energy and momentum are sufficiently 
different to indicate that Sgr and Pal 5 have rather differ- 
ent orbits. Similarly, Figure [9] indicates that the proper 
motion of Pal 5 is inconsistent with that of the Tl stream 
by 1.3/1.9 mas yr _1 (1.7cr/2.3er) for \i a cos<5 and fj,s re- 
spectively. 

Pal 5 has also been shown to have lengthy (~ 20°), 
well-defined tidal tails of its own (Odenkirchen et al. 
2001, 2003; Rockosi et al. 2002; Grillmair & Dionatos 
2006) and these provide another indicator of the direc- 
tion of motion (although sec Odenkirchen et al. 2009). 
A comparison of the direction of these tidal tails (see, 
e.g., Fig. 11 of Rockosi et al. 2002) to the direction 
of the orbital path of Sgr on the sky shows them to be 
nearly perpendicular to one another (i.e., roughly along 
Bq at constant A©). We therefore conclude that Pal 5 
is unlikely to be a member of the Sgr globular cluster 
system. 

5.13. Pal 12 (ID=53, A = 39° J 

Pal 12 has been the subject of vigorous debate in re- 
cent years, and claims of possible association with the 
Sgr streams have been made by Martincz-Delgado et al. 
(2002), Palma et al. (2002), Bellazzini et al. (2003ab), 
Majewski et al. (2004), and on the basis of full space 
velocities by Dinescu et al. (2000). The chemical abun- 
dance of Pal 12 ([Fe/H] = -0.83; Marin-Franch et al. 
2009, see also Cohen et al. 2004, Armandroff & Da Costa 
1991) and low [a/Fe] signatures characteristic of the Sgr 
family have also led Cohen et al. (2004) and Sbordone 
et al. (2005) to conclude that Pal 12 is likely to be asso- 
ciated with the Sgr stream. 

We find that Pal 12 presents an unusual case, in that it 
has a relatively convincing probability of association with 
both the LI (F 3 ,liy = 0.048) and Tl (P 3 ,tiy = 0.012) 
arms. As illustrated in Figures [6] and this curious cir- 
cumstance arises because Pal 12 lies in a region of the 
South Galactic Cap at which the Tl arm overlaps in 
angular position, distance, and radial velocity with the 
LI arm which has passed through the northern Galactic 
hemisphere, through the Galactic disk, and is complet- 
ing its wrap through the southern Galactic hemisphere. 
The angular position of Pal 12 matches the LI arm best 
(Figure[7]), lying only 5.5° (~ 1.4c) from the centroid of 
the stream, compared to 9° (~ 2.6a) from the centroid 
of the Tl arm (Figures [6] and [8]) . The distance of Pal 12 
(19.1 kpc) best matches the Tl arm though, lying 4 kpc 
(~ 1.4ct) closer than the Tl arm and 11 kpc (~ 3tr) closer 
than the LI arm. The radial velocity is approximately 
equally matched to both segments of the Sgr stream, 
matching the LI arm to within 44 km s _1 (~ 1.3<r) and 
the Tl arm to within 16 km s _1 (~ 1.5<r). 

This situation is not clarified by the observed proper 
motions: As illustrated in Figures 191 and ITUI Pal 12 is dis- 
crepant from the Tl arm by 1.3/1.5 mas yr _1 (3cr/3.5cr) 
in (i a cosS/ns respectively, and is discrepant from the LI 



arm by 0.5/2.7 mas yr _1 (1.4cr/6.2cr). The uncertainties 
given on the proper motion of Pal 12 are particularly 
small (0.3 mas yr _1 ) however; if these have been under- 
estimated the significance of the proper motion discrep- 
ancies decreases commensurately. 

If one supposes that Pal 12 is associated with the Tl 
arm (as favored by the distance and proper motion es- 
timates), it must be explained why the angular position 
of Pal 12 places it at the extreme outer edge of the trail- 
ing stream and why the proper motions deviate from the 
model by a total of ~ 4<r. In contrast, if one supposes 
that Pal 12 is associated with the LI arm (as favored 
by the angular coordinates, as well as by Dinescu et al. 
2000 and Majewski et al. 2004) it must somehow be 
explained why Pal 12 lies at a smaller distance than pre- 
dicted for the LI stream (20 kpc versus 30 kpc) and why 
the proper motions are inconsistent by over 6a. Since 
the LM10 model of the leading tidal stream wrapped so 
far (nearly 360°) from the Sgr dwarf is significantly less 
well constrained than the trailing arm stream in the same 
region of the sky (see, e.g., Majewski et al. 2003), it is 
perhaps more plausible that Pal 12 is associated with the 
LI arm. 

Given the low probability (P3 = 0.012, correspond- 
ing to ~ 0.6 clusters out of a randomly drawn sample 
of 51) for a globular cluster to appear to match the Sgr 
stream as well as Pal 12 by chance, in combination with 
the unique chemical fingerprint characteristic of the Sgr 
globular clusters, it seems likely that Pal 12 is physically 
associated with the Sgr stream. In the future, an in- 
dependent proper motion study of Pal 12 will hopefully 
check which arm of the Sgr stream (if either) Pal 12 is 
associated with. 

5.14. Pal 13 (ID=54, A Q = 70°J 

Pal 13, located towards the South Galactic Cap, has 
a distance and radial velocity consistent with the LI 
wrap of the leading tidal stream (see Figure [7]) • At 15.1 
kpc from the Sgr midplane however (corresponding to 
Bq — —33.5°) it is a poor match to the predicted loca- 
tion of LI tidal debris. The L1Y arm has a mean value 
of Bq = +0.3° with an angular width as,Sgr — 3.2°, 
placing Pal 13 over 10c away from the stream and giving 
a P3 = 37% probability that a randomly chosen artifi- 
cial globular cluster would match the Sgr stream better 
than Pal 13. If the older (> 3 Gyr old) debris from the 
LM10 is included in the analysis, the angular width of 
the stream increases slightly, placing Pal 13 only 6.7cr 
away from the LI stream (i"V,Li = 0.155), but neither 
value is compelling evidence for association. 

Strong evidence against the association of Pal 13 with 
the Sgr stream is provided by the proper motions mea- 
sured by Siegel et al. (2001; Figure ITU]) . which are incon- 
sistent with the n a cosS and fj,$ of the LI arm by 3.2/3.1 
mas yr _1 (7.3<t/4.8it) respectively. In contrast to Bel- 
lazzini et al. (2003a), we therefore concur with Palma 
et al. (2002) that Pal 13 is not associated with the Sgr 
dwarf. 

5.15. Pal 15 (ID=37, A Q = 324° j 

Pal 15 lies on the opposite side of the Galactic Cen- 
ter towards (l,b) = (20°, 25°) at a distance of 44.6 kpc. 
As indicated by Table [3] it has a formal probability 
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P3 = 0.128 for association with the LI arm of the Sgr 
tidal stream. This association is relatively weak given 
the corresponding expectation that ~ 7 globular clusters 
of our sample of 51 would match the Sgr stream better 
than Pal 15 by chance. While Pal 15 is a reasonable 
match to both the distance and radial velocity of the LI 
arm at A Q = 324°, it lies ~ 20° (~ 16 kpc) away from 
the bulk of the LI tidal debris at this orbital longitude. 
Although minor discrepancies in angular position might 
be rationalized in regions where the LM10 model of the 
Sgr streams is poorly constrained (e.g., the L2/T2 arms), 
Pal 15 lies within 35° of the Sgr core and its angular po- 
sition disagrees strongly with observations of the leading 
arm as traced in 2MASS (e.g., Figure 3 of Majewski et 
al. 2003). The strong 2MASS M-giant constraints on the 
width of the Sgr debris stream at this orbital longitude 
rule out Pal 15 as a Sgr globular cluster. 

5.16. Whiting 1 (ID =3, A Q = 103° j 

Whiting 1 (Wh 1) is a relatively poorly known, re- 
cently discovered (Whiting et al. 2002) globular clus- 
ter that is fairly young (~ 6.5 Gyr) and metal rich 
([Fe/H] = -0.65; Carraro 2005). On the basis of a dis- 
tance estimate d = 29.4 kpc and radial velocity measure- 
ment i>gsr = —105 ± 1.8 km s _1 , Carraro et al. (2007) 
suggested that it may be associated with the trailing arm 
of the Sgr stream in the direction of the South Galactic 
Cap. We confirm this association, finding that Whit- 
ing 1 lies within ~ 0.2 kpc of the Sgr plane (closer than 
any other cluster save M 54), and that its angular po- 
sition, heliocentric distance, and radial velocity match 
those of the Tl stream to extremely high accuracy (see 
Figs. Hand HI). We calculate a probability P 3 ~ 0.1% 
(corresponding to 0.05 clusters out of a sample of 51) 
that a globular cluster would match the Sgr stream as 
well as Whiting 1 by chance. While proper motions for 
this cluster are not available to confirm the hypothesis, 
it appears that Whiting 1 can be associated with the Sgr 
system with a similar confidence as for M 54, Arp 2, and 
Terzan 8. 

As indicated by Table [31 there is also a lower confidence 
association of Whiting 1 with the secondary wrap of the 
Sgr leading arm (L2). However, since this association is 
at much lower confidence (P = 6.4%, corresponding to 
an expectation of ~ 3 clusters by chance in our sample 
of 51) than the Tl association (and indeed it is not cer- 
tain if L2 debris is even present in the Galactic halo) we 
therefore believe that the apparent L2 association is not 
significant. 

5.17. M 2, NGC 5466, NGC 5824, NGC 6426, Rup 106 

Bcllazzini et al. (2003a) discuss the possible associa- 
tion of these globular clusters with the Sgr stream, and 
Sbordone et al. (2005) note the chemical similarity be- 
tween Rup 106 and the Sgr family. However, none of 
these clusters present significant evidence for association 
with any wrap of the Sgr arms, with probabilities of 
P 3 = 33%, 53%, 29%, 80%, and 62% that a randomly 
chosen globular cluster in the Galactic halo would by 
chance appear to match the Sgr stream better than M 2, 
NGC 5466, NGC 5824, NGC 6426, or Rup 106 respec- 
tively. While the poor match of the angular positions, 
distances, and radial velocities of these clusters is suffi- 



cient to exclude them as Sgr members, the proper motion 
reinforces this conclusion for M 2. 

5.18. AM 4, BE 176, ESO 280-SC06, IC4499 

Radial velocities have not been measured for these four 
clusters, and their association statistics have therefore 
been estimated by repeating the analysis of Section 2] us- 
ing distance and angular position alone. These statistics 
P3 > 0.36 and P7 > 0.45 for all four clusters, indicat- 
ing that neither are likely to be associated with the Sgr 
stream. As illustrated in Figures|5]and[7](see ID numbers 
25, 29, 32, and 41), ESO 280-SC06 is barely consistent 
with the angular position of the T2 stream, but lies much 
too close to the sun (21.7 kpc vs. 79.2 kpc). Neither BH 
176 nor IC 4499 lie within ~ 25° of any simulated stream 
debris. Carraro (2009) make the case that AM 4 might 
be associated with Sgr based on their newly derived dis- 
tance d — 33^4 kpc. However, the closest angular match 
to AM 4 is provided by the LI arm (whose centroid lies 
24.6° away, or 3.7<r), and both the value d = 29.9 kpc 
quoted in Table [5] and the revised estimate of Carraro 
(2009) are strongly inconsistent with the distance of 53 
kpc expected for the LI arm at this angular position. 
While radial velocities for these four clusters would ver- 
ify our analysis, we strongly doubt that any of them is 
associated with the Sgr stream. 

5.19. Globular Clusters at Other Galactocentric Radii 

As described in [J3]above, we have focused our attention 
on globular clusters in the range of Galactocentric radii 
8 < roc < 100 kpc as a generous interval surrounding 
the predicted locations of Sgr tidal debris in the LM10 
model. However, it is conceivable that clusters currently 
located at other radii may also have originated in the Sgr 
dwarf and have been missed in our analysis due to the 
effects of dynamical friction or other limitations of the 
LM10 model. 

While dynamical friction is not expected to substan- 
tially alter the path of individual globular clusters in the 
halo of the Milky Way (the inspiral time for a cluster at 
a radius of 8 kpc in the Galactic potential is expected to 
be much longer than a Hubble time; see, e.g., Eqn 8-24 
of Binney & Tremaine [2008, p. 655]), it is expected to 
have an appreciable effect on the path of the Sgr dSph 
itself which was not incorporated in the LM10 model. 
Since Sgr should have sunk deeper into the gravitational 
potential of the Milky Way over its lifetime, tidal de- 
bris (and any corresponding globular clusters) torn from 
Sgr at early times in its interaction might be expected 
to lie at systematically larger distances than more recent 
epochs of debris. Only two globular clusters are found in 
the Harris (1996) catalog at distances tqc > 100 kpc; E 
1 and Pal 4 at distances of tqc = 123 and 112 kpc respec- 
tively. While El lies ~ 100 kpc from the Sgr plane, Pal 
4 is only ~ 15 kpc away from this plane. At A© = 239°, 
it is perhaps plausible that Pal 4 may be linked to the 
Sgr stream in a similar manner to NGC 2419 (see also 
Palma et al. 2002). If the Tl arm continues to rise in dis- 
tance for Aq > 160° instead of turning around as seen in 
the LM10 model, this arm may pass through both NGC 
2419 and Pal 4. This picture is perhaps unlikely, how- 
ever, given the recent confirmation of the Tl stream at 
A Q ~ 225° with properties similar to those of the LM10 
model predictions (see discussion in £|6.5[) . 
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Since the Milky Way gravitational potential adopted 
by LM10 is imperfect, it is also possible that some 
Sgr tidal debris may fall inside tqc = 8 kpc. Indeed, 
three globular clusters at such radii (M 5, NGC 6356, 
and Terzan 3) have been discussed by previous authors 
(Palma et al. 2002; Bellazzini et al. 2003a) as poten- 
tial Sgr cluster candidates based either on their near- 
alignment with the Sgr plane or on their similarity of 
orbital pole families. NGC 6356 in particular lies on the 
edge of our selection annulus (roc = 7.6 kpc) at an- 
gular coordinates (A ,B @ = (336.6°, -5.2°)). With a 
heliocentric distance of d = 15.2 kpc it is most closely 
matched by the L2 wrap of the leading arm, which is 
expected to lie at d w 25 kpc. However, the radial ve- 
locity of NGC 6356 is i>gsr = +64 km s _1 , which is 
inconsistent with the LM10 prediction of the L2 arm by 
more than 250 km s _1 . If we redo our construction of a 
statistical comparison sample from [J3] - 14] to incorporate 
globular clusters at tqc < 8 kpc, we find that none of 
the clusters in this radius range have association statis- 
tics better than P3 or P-j ~ 30%. To within the limi- 
tations of the LM10 model, we are therefore confident 
in our decision to exclude these clusters from detailed 
consideration. 

6. ANALYSIS OF INDIVIDUAL ULTRAFAINT SYSTEMS 

In Figures [6] and [7] we overplot each of the nine ultra- 
faint Milky Way satellites from Table [2] against the LM10 
model of angular positions, distances, and radial veloci- 
ties for the Sgr trailing/leading arm streams respectively. 
As indicated by Figures 0] and [5j the ultrafaint satellites 
are more likely to appear to match the Sgr stream by 
chance than the globular cluster sample by virtue of the 
footprint of the sky area that has been surveyed by the 
SDSS. Although Figure [5] suggests that there is no popu- 
lation of these satellites obviously associated with the Sgr 
stream, for consistency with the globular cluster analy- 
sis we have listed satellites with P < 0.15 in Table [3] 
and P > 0.15 in Table |H and discuss all nine satellites 
individually below. 

6.1. Bodtes I (ID =62, A Q = 277° ) 

Bootes I (Boo I; Belokurov et al. 2006b) is located 
at (l,b) = (358.1°, 69.6°) at a distance of 62 kpc with a 
radial velocity i>gsr = 103 ± 3.4 km s" 1 (Munoz et al. 
2006; Martin et al. 2007). As illustrated in Figures M 
and [71 there is no particularly remarkable match of Boo 
I to any segment of the Sgr stream, with P3 = 0.493 
and P7 = 0.774. While Boo I lies in projection against 
the T1/T2 arms, the distance is mismatched by 41/19 
kpc respectively, and the radial velocity is mismatched 
by 144/262 km s _1 . Similarly, Boo I lies near the edge 
of the L1/L2 streams but is mismatched in distance by 
11/51 kpc and in radial velocity by 149/266 km s _1 . We 
therefore conclude that Boo I is not associated with the 
Sgr stream. 

6.2. Bodtes II (ID=61, A Q = 277°/ 

Bootes II (Boo II; Walsh et al. 2007) is located in the 
direction (l,b) = (353.8°, +68.9°), and has a distance of 
46 kpc and radial velocity vgsr = —116 km s _1 (Koch 
et al. 2009). As illustrated in Figure [3 Boo II matches 
the LI arm to within ~ 11° in Bq, ~ 70 km s _1 in ra- 
dial velocity, and ~ 5 kpc in distance. The Sgr stream is 



expected to be relatively broad at this location however 
(i.e., near the apogalacticon of its leading orbit), so these 
offsets correspond to 2a, 3cr, and la respectively. While 
Boo II is therefore moderately well matched to the Sgr 
stream, the statistical significance of the match is under- 
whelming with P 3 = 0.124 and P 7 = 0.239. That is, in 
a sample of 9 ultra-faint satellites drawn from a random 
distribution about the Milky Way, roughly one to two 
such satellites (depending on whether the P3 or P7 crite- 
ria are used) would be expected to match the properties 
of the Sgr stream as well as Boo II by chance. We there- 
fore concur with Koch et al. (2009) that Boo II is perhaps 
the most likely of the ultra-faint dwarf satellites of the 
Milky Way discovered to date to be associated with Sgr, 
but caution that the case for association is statistically 
very weak. 

6.3. Coma Berenices (ID=60, A Q = 252°J 

Coma Berenices (Coma Ber) was discovered by Be- 
lokurov et al. (2007) at a distance of 44 ± 4 kpc in the 
direction (l,b) = (241.9°, +83.6°), and has a radial ve- 
locity of w GS r = +81.8 ± 0.9 km s" 1 (Simon & Geha 
2007). Belokurov et al. (2009) note that Coma Ber is 
superimposed on the edge of the Sgr stream at a dis- 
tance that suggests association with the 'C wrap of Sgr 
debris in the model of Fcllhauer et al. (2006). This 'C 
stream corresponds to the L2 wrap of the LM10 model; 
as shown in Figure [3 no such association is seen with 
the LM10 model with P 3 = 0.446 and P 7 = 0.434. While 
Coma Ber lies in the midst of the T1/T2 wraps and at 
the outer edge of the L1/L2 wraps, it is mismatched to 
the T1/T2/L1/L2 wraps respectively by 23/14/3/34 kpc 
in distance and 26/225/182/136 km s _1 m radial veloc- 
ity. We therefore conclude that Coma Ber is probably 
not associated with the Sgr system. 

6.4. Pisces I (ID=64, A Q = 73°; 

Pisces I was originally identified as RR Lyrae 'Struc- 
ture J' in the SDSS Stripe 82 catalog by Sesar et al. 
(2007), and recently spectroscopically confirmed as an 
ultra- faint dwarf galaxy by Kollmeier et al. (2009). Ly- 
ing in the direction (l,b) = (88°,— 58°), Pisces I lies at 
the outer edge of the L2 debris stream, but its large dis- 
tance (d = 85 kpc) is inconsistent with identification with 
anything but the T2 wrap of the Sgr stream. While the 
radial velocity of Pisces I («gsr = 45.4 ± 4 km s _1 ) is 
consistent with the T2 identification, its angular position 
is discrepant with the LM10 model of the T2 wrap by 
~ 37° (~ 4a). Since Pisces I therefore has formal associ- 
ation statistics P3 = 0.999 and P7 = 0.290, we conclude 
that it is unlikely to be associated with Sgr. 

6.5. Segue 1 (ID=58, A Q = 225°/ 

Segue 1 was discovered in the SDSS by Belokurov et 
al. (2007), and lies at a distance of 23 ± 2 kpc in the 
direction (I, b) = (151.763°, +16.074°). The radial veloc- 
ity of Segue 1 is v G sr = 114 km s" 1 (Geha et al. 2009). 
The nature of Segue 1 is a matter of debate: While Be- 
lokurov et al. (2007) favor the interpretation that it is 
an unusually extended globular cluster, the velocity dis- 
persion measured by Geha et al. (2009; a = 4.3 ± 1.2 
km s _1 ) leads these authors to conclude that Segue 1 
is a dark-matter dominated dwarf galaxy. Contradictory 
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claims have also been made both for (Belokurov et al. 
2007; Niederste-Ostholt et al. 2009) and against (Geha 
et al. 2009) the association of Segue 1 with the Sgr tidal 
tails. 

If Segue 1 is associated with the Sgr stream, Table 
[3] suggests that the most likely such association is with 
the Tl trailing wrap of debris wrapped into the northern 
Galactic hemisphere (the radial velocity is mismatched 
with the T2/L1 arms by 190/249 km s _1 respectively, 
and the distance is a factor of ~ 2 too great to match 
the L2 arm). While the distance and radial velocity of 
Segue 1 are both relatively good matches to the Tl wrap 
(see Figure |6]), the angular position of Segue 1 differs 
from the centroid of the stream by 18.6°. If only the 
most recent 3 Gyr of tidal debris (which best correspond 
to epochs for which there is conclusive observational evi- 
dence) are considered, we should expect the T1Y stream 
to be relatively narrow at this point (crs.Sgr = 4.3°), so 
that Segue 1 is offset by ~ 4.3<r from the T1Y stream. 
If tidal debris from > 3 Gyr ago is included, the angular 
width of the stream increases to cs.Sgr — 8°, decreasing 
the discrepancy to only ~ 2.3cr. Formally, we calculate 
Segue 1 to have association statistics P3 = 0.318 and 
P7 = 0.113. Neither is particularly compelling evidence 
for association with the Sgr stream since ~ 3 such satel- 
lites (using the P3 statistic) in a given sample of nine 
drawn from a randomly distributed halo population are 
expected to match the stream better than Segue 1 simply 
by chance. 

We note, however, that Niederste-Ostholt et al. (2009) 
may offer an important confirmation of the LM10 model 
of the Sgr stream at the orbital longitude of Segue 1 
(Aq = 225°), bolstering our confidence in the LM10 pre- 
diction for the location and properties of the wrapped Tl 
trailing arm. In particular, their Figure 10 plots t>GSR vs. 
declination for blue horizontal branch (BHB) stars from 
the SDSS in the range 145° < a < 155° and shows two 
velocity peaks at «gsr = —105 km s _1 and vgsr = 130 
km s _1 . The LM10 model suggests that these features 
may both be due to Sgr tidal debris, the negative velocity 
feature matches that expected for the LI arm to within 
20 km s" 1 , and the positive velocity feature matches that 
expected for the Tl arm to within 2 km s^ 1 . This agree- 
ment is particularly notable for the Tl arm since this may 
demonstrate that the LM10 model is consistent with ob- 
servations at least 60° further along this arm than has 
been previously constrained by the M-giant observations 
of Majewski et al. (2004). However, it is also possible 
that the positive-velocity feature may instead trace the 
Orphan stream (see, e.g., Newberg et al. 2010) which 
is expected to have a similar velocity signature in this 
region of the sky. 

Although the Tl arm may be visible in Figure 10 of 
Niederste-Ostholt et al. (2009), it (as indicated by the 
LM10 model) is centered around S « 25°. Within a 
small field (< 1°) centered on Segue 1 at 8 = 16°, no 
Tl stars should be visible and the velocity signature of 
Sgr in this small field is expected to be entirely that of 
the LI arm (wqsr = —105 km s _1 ). In contrast to the 
suggestion of Niederste-Ostholt et al. (2009), it therefore 
seems unlikely that recent measurements of the velocity 
dispersion of Segue 1 (Geha et al. 2009) have been biased 
by Sgr stream stars. 



Some claims have also been made for the detection 
of tidal tails emanating from Segue 1 (Belokurov et al. 
2007; although c.f. Martin et al. 2008) oriented almost 
perpendicular to the Sgr stream. If these putative tails 
lie along the direction of motion, their orientation may 
further strengthen the conclusion that Segue 1 is not as- 
sociated with Sgr. However, at present the tails extend to 
such small angular distances (~ 2r g; Geha et al. 2009) 
that it is not obvious that they are necessarily genuine 
or good indicators of the orbital path of the dwarf (see, 
e.g., Johnston et al. 2002; Montuori et al. 2007; Klimen- 
towski et al. 2009; Odenkirchen et al. 2009). Conclusive 
confirmation or rejection of the hypothesis that Segue 1 
is associated with the Sgr tidal stream may await the 
accurate determination of Segue l's proper motion and 
comparison to the values \i a cos5 = —2.4 ±0.4 mas yr^ 1 
and [ig = — 1.3±0.1 mas yr _1 expected for the Tl stream 
at this Aq. At present, however, the case for Sgr asso- 
ciation of Segue 1 does not appear to be statistically 
compelling. 

6.6. Segue 2 (ID=56, A = 119°; 

Segue 2 (Belokurov et al. 2009) lies at a distance of 35 
kpc in the direction (/, b) = (149.4°, -38.1°) with a radial 
velocity of i>gsr — 43.3 km s _1 . The angular position, 
distance, and radial velocity of Segue 2 are overplotted 
on the Sgr leading and trailing tidal streams in Figures 
[6] and we concur with Belokurov et al. (2009) in not- 
ing that Segue 2 lies near the edge of the Tl Sgr stream 
at (A ,B ) = (119°, -17°). While Segue 2 also lies at 
roughly the same distance as the Tl stream, its radial 
velocity is inconsistent with this identification by 182 
km s _1 (17cr). Indeed, we calculate a P3 = 89.2% and 
Pi = 33.0% chance that an artificial ultra-faint satel- 
lite drawn from a random angular distribution about the 
Milky Way would match the Sgr streams better than 
Segue 2 by chance, leading us to strongly doubt that 
Segue 2 is associated with the Sgr stream. 

6.7. Segue 3 (ID=63, A = 41°; 

Segue 3 (Belokurov et al. 2010) has recently been iden- 
tified as an ultra- faint star cluster at a distance of 16 kpc 
in the direction (l,b) = (69.4°, -21.3°). Although no 
radial velocity measurement has yet been made of this 
system, we estimate an association statistic P3 = 0.945 
based on its distance and angular position. This ex- 
tremely high value reflects the fact that Segue 3 lies > 30° 
away from any segment of the Sgr stream (indeed, at 
Bq — —48° it does not even fall within the plot window 
of the uppermost panel of Figures [6] and [7)). We therefore 
strongly doubt that Segue 3 is associated with Sgr. 

6.8. Ursa Major II (ID =57, A = 202°; 

Ursa Major II (UMa II; Zucker et al. 2006) lies at 
(l,b) = (152.5°, +37.4°) with a distance of 32 kpc and 
a radial velocity «gsr = —33.4 ± 1.9 km s _1 (Simon & 
Geha 2007). As illustrated by Figure El UMa II lies at 
(A Q ,B e ) = (202°, -33°) (about 3er away from the T2 
stream at Bq = —10.9° ± 7.0°), and is roughly consis- 
tent with the distance (25.2±1.4 kpc) and radial velocity 
(—8.6 ± 16 km s" 1 ) of this segment of the stream, es- 
pecially considering the large uncertainties in the LM10 
model wrapped this far from the Sgr dSph. However, this 
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rough consistency is not significant relative to the typical 
apparent correlation between artificially generated ultra- 
faint satellites and the Sgr streams, with P3 = 0.954 and 
P-i = 0.185. We are therefore inclined to doubt the asso- 
ciation between UMa II and Sgr, although future proper 
motion measurements and observational constraints on 
the Sgr stream at such large angular separations from 
the dSph will be required in order to conclusively decide 
the matter. If the LM10 model is still accurate for the 
T2 arm wrapped ~ 360° beyond the last conclusive ob- 
servational constraint, we might expect the T2 stream 
at the location of UMa II to have a proper motion \x a 
cos^ = —2.2 ± 0.2 mas yr _1 and us — —2.4 ± 0.2 mas 
yr _1 . 

6.9. Willman I (ID=59, A Q = 223°J 

Willman I (Willman et al. 2005) lies at (l,b) = 
(158.57°, +56.78°) with a distance of 38 kpc and a ra- 
dial velocity wqsr = 35.4 ± 2.5 km s^ 1 (Martin et al. 
2007). As for Boo I and Coma Ber, Willman I is not ob- 
viously associated with any wrap of the Sgr stream with 
P 3 = 0.940 and P 7 = 0.694. In Figures M and we over- 
plot Willman I on the LM10 model of the Sgr stream, 
noting that the satellite lies ~ 20° away from the L1/L2 
wraps. While Willman I lies nearer to the T1/T2 wraps, 
it is more distant than anticipated for the T1/T2 streams 
(28/26 kpc respectively) and has a radial velocity incon- 
sistent with these wraps by 102/101 km s _1 respectively. 
We therefore conclude that Willman I is probably not 
associated with the Sgr system. 

6.10. Satellite Systems at Large Galactocentric 
Distances 

We have omitted from consideration above the ultra- 
faint satellites Canes Venatici I, Canes Venatici II, Her- 
cules, Leo IV, Leo V, Leo T, Pisces II, and Ursa Ma- 
jor I since these systems lie at galactocentric distances 
tqc > 100 kpc at which the LM10 models predicts no 
tidal debris. However, as discussed previously in £ 15. 191 
the LM10 model does not take into account the effects of 
dynamical friction on the orbit of Sgr, and it is possible 
that stars or stellar subsystems lost from the satellite at 
early times may lie at significantly larger distances than 
predicted by the LM10 model. Ultra-faint dwarf satel- 
lites may also not have actually formed in the potential 
well of the Sgr dwarf, but may have been gravitationally 
associated with it and fallen into the Milky Way as a 
group (e.g., Li & Helmi 2008; D'Onghia & Lake 2009). 
Such associated systems might be expected to be among 
the first casualties of tidal stripping by the Milky Way 
and therefore lie at large distances relative to the present- 
day Sgr core and recent tidal debris. 

All of these satellites lie out of the present orbital 
plane of Sgr however: Canes Venatici I (|^Sgr| = 86 
kpc; Bq = —23°), Canes Venatici II (|^Sgr| = 53 kpc; 
B Q = 20°), Hercules (|Z Sgr | = 66 kpc, B Q = -28°), Leo 
IV (|Z Sgr | = 50 kpc, B Q = 19°), Leo V (|Z Sgr | = 53 
kpc, B Q = 17°), Leo T (|Z Sgr | = 84 kpc, B Q = 12°), 
Pisces II (|Zsgr| = 86 kpc, B Q = 28°), and Ursa Ma- 
jor I (|Z Sgr | = 46 kpc, Bq = 25°) all have |Z Sgr | > 45 
kpc and Bq > 10°. Similarly, of the other Local Group 
dwarf galaxies only Leo T, Leo A, Leo I, Leo II, and the 
Sagittarius dwarf irregular galaxy (Sgr DIG) lie within 



TABLE 5 
Star Clusters Associated with the Sgr Stream 



ID 


Name 


Arm 


High Confidence 


3 

27 
12 
45 
46 


Whiting 1 
NGC 5634 
M 54 
Arp 2 
Terzan 8 


Tl 
Tl 
CORE 
Tl 
Tl 


Moderate Confidence 


11 

22 
11 
53 


Berkeley 29 
NGC 5053 
Terzan 7 
Pal 12 


LI 

Tl 

Tl 

Tl/Ll 


Low Confidence 


7 
18 


Pal 2 

NGC 4147 


Tl 
Tl 



15° of the Sgr orbital plane; all of these systems lie be- 
tween 200 kpc and 1 Mpc distant from the Milky Way. 
If substantial evolution has occurred in the orbit of the 
Sgr dSph over its lifetime, it is conceivable that some of 
these systems may eventually prove to be associated in 
some manner with Sgr (whether forming inside the dSph 
or gravitationally associated with it). However, it is not 
possible to meaningfully constrain this eventuality at the 
present time using the LM10 model. 



7. 



DISCUSSION: THE SGR CONTRIBUTION TO THE 
GALACTIC HALO 



In Table [5] we summarize the star clusters that we con- 
clude may be associated with Sgr to high/moderate/low 
confidence, and the tidal arm in which they appear to re- 
side. Out of 15 candidate Sgr globular clusters selected 
via the criterion P 3 < 0.15, we find that five (Arp 2, 
M 54, NGC 5634, Terzan 8, and Whiting 1) are very 
likely to be associated with the Sgr dwarf, an additional 
4 (Berkeley 29, NGC 5053, Pal 12, Terzan 7) are moder- 
ately likely to be associated, two (NGC 4147, Pal 2) may 
be associated with relatively low confidence, and four (M 
53, NGC 288, Pal 5, Pal 15) are very unlikely to be as- 
sociated. Taking as our final Sgr globular cluster sample 
those in the high- and moderate-confidence categories, 
we conclude that 5-9 globular clusters are genuinely as- 
sociated with Sgr while 6-10 simply appear to be asso- 
ciated by chance. This conclusion is consistent with our 
expectation, based on statistical realizations of a ran- 
domly distributed artificial globular cluster population, 
that 8 ± 2 of the 15 candidates would be genuine Sgr 
clusters. As illustrated by Table [TJ the clusters which 
we conclude are most likely to be associated with Sgr 
are generally those for which there was the most consis- 
tent agreement on Sgr membership in the literature. We 
are able, however, to rule out various candidates such 
as M 2, NGC 288, and NGC 6356 which have also been 
proposed by various authors. Our overall number of Sgr- 
member globular clusters is consistent with Bellazzini et 
al. (2003a), who concluded (at 95% confidence) that at 
least eight globular clusters are likely to be associated 
with the Sgr system; four 'outer-halo' globular clusters 
in addition to the four 'core' clusters M 54, Arp 2, Terzan 
7, and Terzan 8. 

Curiously, of the 5-9 globular clusters which we as- 
sociate with the Sgr stream only one or two (Berkeley 
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Fig. 11. — Vital diagram (after Gnedin & Ostriker 1997) for 
high, moderate, and low- confidence Sgr globular clusters. Masses 
and half-light radii for the classical globular clusters are taken from 
tabulations by Gnedin & Ostriker (1997) and Harris (1996). These 
parameters are estimated for Berkeley 29 and Whiting 1 using data 
from Kaluzny (1994), Tosietal. (2004), and Carraro (2005). Over- 
laid in black solid lines are the survival regions for distances roc = 
12, 8, and 5 kpc from the Galactic center in the Ostriker & Caldwell 
(1983) galaxy model. 

29, and possibly Pal 12) lie in the leading arm of the 
Sgr stream. In part, this may be due to selection bias; 
~ 75% of the clusters assigned to the Tl arm lie within 
~ 100° of the Sgr dwarf. The corresponding section of 
the LI arm lies almost entirely in the zone of avoidance 
at low Galactic latitudes where sparse clusters such as 
Whiting 1 and Arp 2 would be more difficult to detect. 
However, this may also point to a difference in the sur- 
vivability of clusters in the leading versus trailing Sgr 
tidal arms. Since clusters in the leading arm would pass 
more closely to the Galactic center than clusters in the 
trailing arm, gravitational shocks from the Galactic disk 
may be more efficient in disrupting clusters in the leading 
arm. Indeed, as illustrated in Figure [Til sparse clusters 
such as Whiting 1 in particular lie significantly outside 
the survival regime for orbits 12 kpc from the Galactic 
center (roughly the perigalacticon of the Sgr leading arm 
tidal debris on some orbits in the triaxial halo model of 
LM10), and may therefore be disrupted relatively quickly 
once they are removed from Sgr into the leading tidal 
stream. 

It is perhaps unsurprising that we are unable to asso- 
ciate any stellar clusters with the older T2/L2 wraps of 
Sgr tidal debris confidently, given the possible destruc- 
tion of these clusters over time, the current lack of evi- 
dence for these sections of Sgr stream in the Galactic halo 
and the correspondingly poor constraints on the LM10 
model in this regime. Our null result may, however, con- 
stitute further evidence that the Sgr dwarf has only been 
losing appreciable quantities of stellar debris in the last 
~ 3 Gyr, which points to a constraint on the interac- 
tion time of the dwarf similar to that currently implied 
by wide-field observations of the stellar debris streams. 
Alternatively, discrete stellar systems such as globular 
clusters may simply have formed sufficiently deep in the 
Sgr gravitational potential well that they have only been 



stripped from the dwarf after it has already experienced 
significant tidal disruption. 

We find no compelling evidence to associate Sgr with 
any of the population of ultra-faint Milky Way satel- 
lites recently discovered in the SDSS. While some of the 
ultra-faint satellites (e.g., Boo II, Segue 1) are roughly 
co-spatial with sections of the Sgr stream and have sim- 
ilar radial velocities, this apparent correlation may be a 
product of the survey area in which the satellites were 
discovered; the ultra-faint satellite population is statis- 
tically consistent with being unassociated with the Sgr 
tidal streams. It is not possible for us to discriminate 
between the cases in which some of these systems (a) 
formed within Sgr, (b) were loosely associated with Sgr 
and fell into the Milky Way as a group (e.g., Li & Helmi 
2008; D'Onghia & Lake 2009), or (c) are completely un- 
related to Sgr. We are therefore unable to cast any light 
upon the origins of these enigmatic systems, except that 
they were not likely previously tightly bound parts of the 
Sgr dSph. 

Niederste-Ostholt et al. (2010) have recently provided 
an estimate of the absolute magnitude of the initial Sgr 
dwarf by summing the light contained in the remnant 
Sgr core and the tidal streams. Adopting the fainter 
of their two luminosity functions (which is most consis- 
tent with the LM10 assumption that the present abso- 
lute magnitude of the Sgr core My — —13.64), assuming 
that the tidal tails are symmetric in mass content, and 
assuming that there is negligible luminosity contained 
in regions of the stream for which there is no observa- 
tional evidence in the combined 2MASS + SDSS sur- 
veys, the Niederste-Ostholt et al. (2010) analysis in- 
dicates that the original Sgr dwarf had My rj —15.0 
(with ~ 70% of this light now residing in the tidal tails). 
Our 5-9 globular clusters therefore correspond to a spe- 
cific globular cluster frequency Sn — 5 — 9 for the orig- 
inal Sgr dSph, where S N = Nqc x 10 04 ( M ^+ 15 ) (Harris 
& Van den Bergh 1981). While this frequency is high 
compared to that of the Milky Way, LMC, and SMC 
(Sn = 0.7, 0.8, 1.0 respectively based on the compilation 
by Forbes et al. 2000), it is substantially less than that 
of Fornax (Sn = 29) and of previous estimates of the 
specific frequency of Sgr (Sn ~ 25; e.g., Forbes et al. 
2000; Forbes & Bridges 2010)0 and similar to that of 
many dE galaxies (Strader et al. 2006; Brodie & Stradcr 
2006, and references therein). 

We caution, however, that our analysis is subject to 
numerous uncertainties, the greatest of which is the un- 
certainty in the LM10 model of the Sagittarius tidal 
streams. While the LM10 model is reasonably well con- 
strained for the most recent 3 Gyr of tidal debris for 
which there is the strongest evidence to date in the 
Galactic halo, its predictions for the possible locations 
and properties of earlier epochs of tidal debris are only 
extrapolations. If significant evolution has occurred in 
the orbit of Sgr over time (e.g., via dynamical friction) it 
may be difficult to identify debris (whether stars or larger 
stellar systems) from the earliest epochs whose present 
orbits may not resemble the current orbit of Sgr. Our 
inventory of the substructures associated with Sgr may 

12 The main reason for the discrepancy between our estimate 
of the specific frequency of Sgr and that of previous studies is our 
assumption of a significantly brighter original satellite. 
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therefore be an underestimate. In future, proper mo- 
tions of the quality expected to be provided by NASA's 
Space Interferometry Mission (SIM) are expected to per- 
mit us to assess the case for membership of individual 
stellar systems with the Sgr tidal stream with greater 
confidence. 

7.1. The Age-Metallicity Relation and Horizontal 
Branch Morphology of Sgr Globular Clusters 

In Figure [T^] (left panel) , we plot the ages and metal- 
licities (drawn from tabulations by Marin-Franch et 
al. 2009, and Forbes & Bridges 2010) of our eleven 
high/moderate/low confidence Sgr globular clusters. As 
discussed recently by Forbes & Bridges (2010) for a 
sample of Sgr-candidate clusters drawn from the lit- 
erature, the clusters describe a relatively well defined 
age-metallicity relation (AMR), which we note is gen- 
erally consistent with the closed-box Sgr AMR derived 
by Siegel et al. (2007). We identify three groupings of 
clusters in particular in Figure [H M 5^3, NGC 5053, 
NGC 5634, Pal 2, and Ter 8 resemble the oldest Sgr 
stellar population ('M54 MPP'), Arp 2 and NGC 4147 
resemble the second-oldest Sgr stellar population ('Sgr 
MPP'), and Be 29, Pal 12, Ter 7, and Whiting 1 resem- 
ble the intermediate age Sgr stellar population ('SInt'). 

In Figure [12] (right panel) we show the Lee diagram 
of the candidate Sgr globular clusters, plotting metallic- 
ity versus horizontal branch (HB) morphology HBR = 
{B -R)/(B + V + R), where B, V, and R represent the 
number of blue, variable, and red horizontal branch stars 
respectively. As has often been noted previously (e.g., 
Mottini & Wallcrstcin 2008, and references therein), the 
Sgr candidate globular clusters are a bimodal group of 
red- and blue-HB clusters. Intriguingly, the strongest 
Sgr cluster candidates (i.e., cyan filled squares) are those 
which, despite their ~ 6 Gyr spread in ages, show the 
least evidence for a second-parameter affecting their HB 
morphologies and lie closest to the single-parameter (i.e., 
metallicity) locus defined by Zinn (1993a; solid line in the 
right panel of Figure n^jlP^l The moderate-confidence 
clusters Be 29, Pal 12, and Ter 7 also follow the single- 
parameter locus, while NGC 5053 and the two low- 
confidence clusters NGC 4147 and Pal 2 have HBs red- 
der than expected for their metallicity. The absence of a 
second-parameter effect in the Sgr globular clusters is 
in marked contrast to the globular cluster systems of 
the LMC and Fornax (e.g., Zinn 1993b; Buonanno et 
al. 1999; Mackey & Gilmore 2004), which show mod- 
erate/strong influences of a second parameter on their 
respective HB types at a given metallicity. This may sug- 
gest that the strength of the second-parameter effect on 
the HB morphologies of accreted globular clusters could 
somehow be determined by the subhalo in which they 
were formed, offering a criterion by which to associate 
clusters in the halo of the Milky Way with their parent 
satellites. 

Although there is a strong gradient in the chemistry of 
stars along the Sgr streams (see, e.g., Chou et al. 2007, 

13 Note that the age and metallicity of M 54 was adopted from 
the Siegel et al. (2007) M 54 MPP measurement. 

14 While M 54 lies to the left of the single-parameter locus, this 
is likely due to confusion between the blue HB of M 54 and the red 
HB of the more metal-rich Sgr stellar populations. 



2010), there is no comparable gradient for the Sgr glob- 
ular clusters. While the two globular clusters wrapped 
farthest from the Sgr dwarf (NGC 5053 and NGC 5634, 
wrapped 266° and 293° from Sgr respectively) are the 
most metal-poor of the clusters, their closest chemical 
kin (Arp 2, M 54, Ter 8) lie < 10° from Sgr. Similarly, 
the relatively metal-rich clusters Whiting 1 and Berkeley 
29 lie in the dynamically young (103° separation from 
Sgr) and old (183° separation from Sgr) sections of the 
stream respectively. Indeed, if Berkeley 29 formed 4.5 
Gyr ago (Carraro et al. 2004; 3.4 - 3.7 Gyr ago from the 
estimate of Tosi et al. 2004) in the 'SInt' burst of star 
formation, the LM10 model suggests that it must have 
been stripped from Sgr less than a Gyr after its forma- 
tion since it would have taken ^3 — 4 Gyr since stripping 
to reach its present longitudinal separation from Sgr. 

8. SUMMARY 
Our major conclusions may be summarized as follows: 

1. We find that five globular clusters (Arp 2, M 54, 
NGC 5634, Terzan 8, and Whiting 1) are very likely 
to be associated with the Sgr dwarf, an additional 
4 (Berkeley 29, NGC 5053, Pal 12, Terzan 7) are 
moderately likely to be associated, and two (NGC 
4147, Pal 2) may be associated with relatively low 
confidence. Confining our attention to those in the 
high- and moderate-confidence categories, we con- 
clude that 5-9 globular clusters are genuinely as- 
sociated with Sgr, consistent with our expectation 
of 8 ± 2 clusters based on statistical realizations 
of a randomly distributed artificial globular cluster 
population. 

2. Based on updated estimates of the initial luminos- 
ity of the pre-interaction Sgr dwarf (My = —15), 
we estimate a specific globular cluster frequency 
Sn = 5 — 9 typical of dE galaxies. 

3. The globular clusters identified as most likely to be 
associated with Sgr are consistent with the AMR 
of Sgr's stellar populations, and show no evidence 
for a second-parameter effect shaping their hori- 
zontal branch morphologies. If the presence of a 
second-parameter effect is a reliable indicator of 
Sgr membership, it suggests that the two lowest- 
confidence Sgr clusters (NGC 4147 and Pal 2) and 
the moderate-confidence cluster NGC 5053 may 
not be genuinely associated with Sgr. There is 
no obvious correlation of age, metallicity, or hor- 
izontal branch morphology with distance from Sgr 
along the tidal streams. 

4. We find no compelling evidence to associate Sgr 
with any of the population of ultra-faint Milky 
Way satellites recently discovered in the SDSS. The 
ultra-faint satellite population is statistically con- 
sistent with being unassociated with the Sgr tidal 
streams, but we are unable to rule out association 
conclusively in all cases. 

We caution, however, that our results are based upon 
comparisons to sections of the Sgr stream that we 
presently observe in the Galactic halo, and do not ac- 
count for possible variations in the orbit of Sgr over time. 
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Fig. 12. — Left panel: Age-metallicity relation for candidate Sgr globular clusters (where observational data is available). Filled cyan 
squares represent high confidence candidates, green crosses moderate confidence candidates, and red open boxes low confidence candidates. 
Ovcrplotted is the AMR of Siegel et al. (2007; solid line), which best fits their reconstructed star formation history of Sgr (elliptical shaded 
regions). Right panel: Horizontal branch morphology for candidate Sgr globular clusters as a function of metallicity. Overlaid (solid line) 
is the 'first-parameter' locus defined by Zinn (1993a), objects to the left of this line are those with a significant second-parameter effect 
shaping their horizontal branch morphologies. 

understand the accretion origins of substructure within 
our Galactic halo. 



For instance, we do not address here the (likely) possi- 
bility that dynamical friction has altered the orbit of Sgr 
over its lifetime, nor are we able to constrain the pos- 
sibility that some objects may have been loosely affili- 
ated with Sgr and fallen into the Milky Way as a group. 
Future observational constraints on the length and char- 
acteristics of the Sgr tidal streams will help to further- 
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